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The status of thermal model descriptions of particle production in heavy 
ion collisions is presented. We discuss the formulation of statistical mod- 
els with different implementation of the conservation laws and indicate 
their applicability in heavy ion and elementary particle collisions. We 
analyze experimental data on hadronic abundances obtained in ultra- 
relativistic heavy ion collisions, in a very broad energy range starting 
from RHIC/BNL (^/i = 200 A GeV), SPS/CERN {yfs ~ 20 A GeV) 
up to AGS/BNL (Vi ~ 5 A GeV) and SIS/GSI {yfs ~ 2 A GeV) to 
test equilibration of the fireball created in the collision. We argue that 
the statistical approach provides a very satisfactory description of ex- 
perimental data covering this wide energy range. Any deviations of the 
model predictions from the data are indicated. We discuss the unified de- 
scription of particle chemical freeze-out and the excitation functions of 
different particle species. At SPS and RH1C energy the relation of freeze- 
out parameters with the QCD phase boundary is analyzed. Furthermore, 
the application of the extended statistical model to quantitative under- 
standing of open and hidden charm hadron yields is considered. 



2 



Peter Braun- Munzinger, Krzysztof Redlich, Johanna Stachel 



Contents 



ll. Introduction! 
"" ^T^^l^irlfal conditions in A-A collisions and deconfinemend 



2 



2^^^^^stk^^nogroac^^^rM^^anonira^crai^isrrj|_ 



2j2^^^e«na^n^^i^o^^^cJ^^eW^^2n^A^^^c^^Hy[C^Mr^iesJ_ 



2 
4 
8 
9 
13 



2 : 3 J __Cjmrj3ax^ with thermal model predictions! 20 



^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ 



13. Exact Im plementation of the conservation laws in the statistical models! 



21 
22 



|3.1. 


Kinetics of time evolution and eauilibration of charged particles! 


24 




3.1.1. 


Kinetic master eauation for probabilities] 


25 




3.1.2. 


The eauilibrium solution of the general rate eauation! 


30 




3.1.3. 


The master eauation in the presence of the net chared 


32 




3.1.4. 


The kinetic eauation for different particle species! 


34 


|3.2. 


The canonical description of an internal symmetry - protection method 1 35 




3.2.1. 


Canonical models with a non-Abelian symmetry 1 


40 




3.2.2. 


The canonical partition function for Abelian charges! 


43 




3.2.3. 


The eauivalence of the canonical formalism in the grand canonical limitl 



4^[^^C^n^^Hie^v^^ori^ollrslp^i^at^I^en^r^tes|" 



4.2. 



52 
52 
60 



EX 



Par^^k^rjrodurtipri^i^ig^^ner^^^^^oJIisjpnd^ 

^^^^^^^^^^^^^^^^ MTrzatlpnamc^tAi^Xvnamics in p-p collisions! 64 
^^^^^^^^^^^^^^^^^^^ 68 



4^3T J ^Statistical Recombination Modell 
4.3.2. Results! 



71 



4 J ^yj^Cjh<ffl | noniuinProchjg^ Collisions at LHC Energy! 74 



4 ri 3^^^COTc^rsipn^£nJjl^a^^^UcU^^Prochj£Wa 



15. Unifedconditi£nsof_pa£^ 



6^F^^^£k^i£kl^^i^^heir^ne^^^fi£eiideiice|_ 



7. Lifting of the strangeness suppression in heavy ion collisions! 
References) 



75 
75 
79 
81 
92 
98 
101 
101 



1. Introduction 

The ultimate goal of the physics program with ultrarelativistic nucleus- 
nucleus collisions is to study the properties of strongly interacting matter 
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under extreme conditions of high energy density. Quantum Chromodynam- 
ics (QCD) predictff - ^ that strongly interacting matter undergoes a phase 
transition from a state of hadronic constituents to a plasma of deconfined 
quarks and gluons (QGP). By colliding heavy ions at ultrarelativistic en- 
ergies, one expec ts t o create matter under conditions that are sufficient for 
deconfinemen4iH"^4 Thus, of particular relevance is finding experimental 
probes to check whether the produced medium in its early stage was in- 
deed in the QGP phase. Different probes have been studied with the various 
SPS/CERN and RH1C/BNL experiments. The most promising signals of 
deconfinement are related to par ticular pr operties of the trans verse mo- 
mentum spectra of photon^^IEl dileptonJ-^H-^, and hadrons EEDH22I 
The photon rate is studied to probe the temperature evolution from for- 
mation to decoupling of the fireball, implying sensitivity to a high tem- 
perature deconfined phase. The invariant mass distribution of dileptons 
is expected to be modified by in-medium effects related to chiral symme- 
try restoratio n ™al ( Uifl^ mil The modification of charmonium produc- 
tion wa s arg ued to be a consequence of collective effects in the deconfined 
mediuirPESl. 

Hadron multiplicities and their correlations are observables which can 
provide information on the nature, composition, and size of the medium 
from which they are originating. Of particular interest is the extent to 
which the measured particle yields are showing equilibration. The appear- 
ance of the QGP, that is a partonic medium being at (or close to) lo- 
cal thermal equilibrium and its subsequent hadronization during the phase 
transition should in general drive hadronic constituents towards chemical 
equilibriurrP32ECIIII Conseq uently, a high level of chemical saturation, par- 
ticularly for strange particle a ' \ could be related to the deconfined phase 
created at the early stage of heavy ion collisions. 

The level of equilibrium of secon daries in he avy ion collisions was tested 
by analyzing the par ticle abundance d^ l ^ H^ or their momentum spectra 
| 9 | 20 | 2 1 |-| 23 | 37 | 46 |4T| j n first case one establishes the chemical compo- 
sition of the system, while in the second case additional information on 
dynamical evolution and collective flow can be extracted. 

In this review we will discuss the formulation of statistical models and 
their applications to a phcnomcnological description of particle production 
in nucleus-nucleus collisions. We emphasize the importance of conservation 
laws and their different implementations in the statistical approach. We 
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analyze experimental data on hadronic abundances obtained in ultrarel- 
ativistic heavy ion collisions, in a very broad energy range starting from 
RHIC/BNL (Vs = 130 A GeV), SPS/CERN [y/s ~ 20 A GeV) down to 
AGS/BNL (v^~5A GeV) and SIS/GSI ( v / i~2A GeV) to test equili- 
bration. We argue that the statistical approach provides a very satisfactory 
description of experimental results covering this wide energy range. We 
further provide arguments for a unified description of chemical freeze-out 
of hadrons and discuss excitation functions of different particle species. An 
extension of the model for a quantitative understanding of open and hidden 
charm particle yields will be also discussed. 

1.1. Initial conditions in A-A collisions and deconfinement 

In ultrarelativistic heavy ion collisions, the knowledge of the critical energy 
density e c required for deconfinement as well as the equation of state (EoS) 
of strongly interacting matter are of particular importance. The value of e c 
is needed to establish the necessary initial conditions in heavy ion collisions 
to possibly create the QGP, whereas the EoS is required as an input to 
describe the space-time evolution of the collision fireball a . 

Both of these pieces of information can be obtained today from first 
principle calculations by formulating QCD on the lattice and performing 
Monte-Carlo simulations. In Fig. we show the most recent result of 
lattice gauge theory (LGT) for the temperature dependence of energy den- 
sity and pressure. These results have been obtained in LGT for different 
numbers of dynamical fermions. The energy density is seen in Fig. J2j to ex- 
hibit the typical behavior of a system with a phase transition 13 : an abrupt 
change in a very narrow temperature range. The corresponding pressure 
curve shows a smooth change with temperature. In the region below T c the 
basic constituents of QCD, quarks andgluons, are confined within hadrons 
and here the EoS is well parameterized-^ by a hadron resonance gas. Above 



a In Fig. we show a schematic view of the space-time evolution of heavy ion collisions 
in the Bjorken modeEHl 

b In a strictly statistical physics sense a phase transition in two flavour QCD can only 
appear in the limit of massless quarks where it is of second order. In three flavour QCD, 
with (u,d,s) quarks, the phase transition and its order depends on the value of the quark 
masses. In general it can be a first order, second order or cross-over transition. For 
physical quark masses, both the value of the transition temperature and the order of the 
deconfinement phase transition are still not well established. 
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Fig. 1. Schematic space— time view of a heavy ion collision that indicates four basic 
stages in the evolution of the collision fireball: initial overlap region, pre-equilibrium 
partonic system, equilibrated quark-gluon plasma and its subsequent hadronization to a 
hadron gas. 



T c the system appears in the QGP phase where quarks and gluons can travel 
distances that substantially exceed the typical size of hadrons. The most 
recent results of improved perturbative expansion of the thermodynamical 
potential in continuum QCD indicat e ! I I that, at some distance above T c , 
the EoS of QGP can be well described by a gas of massive quasi-particles 
with a temperature depende nt mas s. In the vicinity of T c the relevant de- 
grees of freedom were arguedSl&SI to be described by Polyakov loops. 

Lattice Gauge Theory predicts, in two-flavour QCD, a critical tem- 
perature T c = 173 ± 8MeV and corresponding critical energy density 
e c = 0.6 ± 0.3 GeV/fm 3 for the deconfmement phase transition^! The 
value of e c is surprisingly low and corresponds quantitatively to the energy 
density inside the nucleon. The initial energy density rea ched in heavy ion 
collisions can be estimated within the Bjorken modePSl. From the trans- 
verse energy Et measured in nucleus-nucleus collisions the initial energy 
density e is determined as 
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Fig. 2. The pressure P and energy density e normalized to the temperature to the 
fourth power, versus temperature normalized to its critical value. The calculation^^! 
were performed within LGT for different numbers of flavors. The values of the corre- 
sponding ideal gas results are indicated by the arrows. 



, s 1 1 dE T 

where the initially produced collision fireball is considered as a cylinder 
of length dz = r dy and transverse radius R ~ A 1 / 3 . Inserting for nR 2 
the overlap area of colliding Pb nuclei together with an assumed initial 
time of Tn ^ 1 fm, and using an average transverse energy at midrapidity 
ecPHat the SPS (y/s = 17.3 GeV) to be 400 GeV, one obtains 



mcasurec 



eg^(r ~ 1 fin) ~ 3.5 ± 0.5 GeV/mT\ (2) 

Increasing the collision energy to ^/s = 130 A-GeV for Au-Au at RHIC 
and keeping the same initial thermalization time as at the SPS, would 
increase eo by only 50-60 %. However, atRHIC the thermalization time 
was argued in terms of different 

modelj™ 

to be shorter by a factor of 

3-5. 

In the context of saturation model jlHI^Z H^I th e thermalization time can 
be possibly related with the saturation scalJ^^l The basic concept of the 
saturation models is the conjecture that there is some transverse momentum 
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scale p sa t where the gluon and quark phase space density saturate J ^ O l ^^ hESl 
For an isentropic expansion of the collision fireball, the transverse energy at 
Psat was related in Ref. I|10fl to that measured in nucleus-nucleus collisions 
in the final state. The saturation scale was also used to fix the thermaliza- 
tion time as r eq ~ l/p S at- Taking the value of p sa t predicted in Ref. |10) for 
RHIC energy p sat ~ 1.13 GeV, one gets r eq ~ 0.2 fm and a corresponding- 
energy density e eq ~ 98 GeV/fm 3 . This is a larger value than expected for 
the initial energy density at RHIC in the McLerran-Venugopalan 

modelMI 

where e RHIC ~ 20 GeV/fm 3 , also in agreement with the prediction of 
Ref. ©• 

At SPS energy the saturation model described in Ref. leads to 
e eq S ~ 16 GeV/fm 3 , a much higher value than that obtained from Eq. J3J). 
The estimate of e eq and initial thermalization time strongly depends on the 
value ofosat and the model assumptions. In a " bottom-up " equilibration 
scenario"' the thermalization time in Au-Au collisions at RHIC energy was 
estimated to be as large as 3.2-3.6 fm and the temperature T ~(210-230) 
MeV. Nevertheless, this initial temperature still corresponds to an energy 
density by factor of 2-3 larger then that required for deconfinement. In A-A 
collisions at the LH C the initial energy density of the equilibrated partonic 
medium is expectecP^ESl to be in the range 400 < t^. q c < 1300 GeV/fm 3 . 

The dominant constituents of the partonic medium produced in ultra- 
relativistic heavy ion collisions at LHC, RHIC and even at SPS energy are 
gluons. The energy density of gluons in thermal equilibrium scales with the 
fourth power of the temperature e = gT 4 , where g is related to the number 
of degrees of freedom. For an ideal gluon gas, g = 167r 2 /30; in an inter- 
acting system, the effective number of degrees of freedom g is smaller. The 
results of LGT shown in Fig. [5] indicate deviations from the Boltzmann 
limit by 20-25 %. Relating the thermal energy density with the initial 
energy density discussed above, one can make an estimate of the initial 
temperature reached in heavy ion collisions. For the SPS, RHIC and LHC 
energies this gives a temperature in the range: 200 MeV < T SPS < 330 
MeV, 210 MeV < T RHIC < 600 MeV and 1000 MeV < T LHC < 1200 
MeV, respectively. 

Comparing the initial energy density expected in heavy ion collisions 
with LGT results, it is clear, that the initial energy density at LHC and 
RHIC by far exceeds the critical value. A large energy density is, however, 
still not sufficient to create a QGP. The distribution of initially produced 
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gluons is very far from being thermaLthus the system needs enough time 
to equilibrate. Recently it was showrPI in the framework of perturbative 
QCD and kinetic theory, that the equilibration of partons should definitely 
happen at LHC and most likely at RHIC energy. 

A previous microscopic studjEa within the Parton Cascade Model has 
led to the conclusion that thermalization can also be reached even at the 
lower SPS energy. Here, however, due to the relatively low collision energy, 
it is not clear whether a model inspired by perturbative QCD is indeed 
applicable. 

Assuming QGP formation in the initial state in heavy ion collisions one 
expects that the thermal nature of the partonic medium could be preserved 
during hadronization. c Consequently, the particle yields measured in the 
final state should resemble a thermal equilibrium population. 

2. Statistical approach - general remarks 

In the approach of Gibbs (see, e.g., Ref. (94)) the equilibrium behavior of 
thermodynamical observables can be evaluated as an average over statisti- 
cal ensembles (rather than as a time average for a particular state). The 
equilibrium distribution is thus obtained by an average over all accessible 
phase space. Furthermore, the ensemble corresponding to thermodynamic 
equilibrium is that for which the phase space density is uniform over the 
accessible phase space. In this sense, filling the accessible phase space uni- 
formly is both a necessary and sufficient condition for equilibrium. Conse- 
quently, the agreement between observables and predictions using the sta- 
tistical operator imply equilibrium (to the accuracy with which agreement 
is observed). "Filling phase space" is not a different statement, although it 
is often and erroneously used in the literature. 

In our further analysis we use in the statistical operator as Hamilto- 
nian that leading to the full hadronic mass spectrum. In some sense this 
is synonymous with using the full QCD Hamiltonian. The only param- 
eters in the statistical operator describing the grand-canonical ensemble 
are temperature T and baryon chemical potential /ib- There is no room 
here for strangeness suppression (7^) factors. So the interpretation is that 
agreement between data and theoretical predictions implies statistical equi- 



c The fact that the phase transition is a driving force towards equilibration is founc 
e.g. in different kinetic models for QGP evolution and hadronization. 
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librium at temperature T and chemical potential fiB- If an additional factor 
7 S is needed to describe the data this implies a clear deviation from chem- 
ical equilibrium: a state in which e.g. strangeness is suppressed compared 
to the equilibrium value implies additional dynamics not contained in the 
statistical operator and not consistent with uniform phase space density. 
Similar arguments, of course, apply if one uses canonical phase space. If, 
in this regime, canonically calculated particle ratios agree with those mea- 
sured, this implies equilibrium at temperature T and over the canonical 
volume V. To the extent that this describes data for e+e- or pp collisions, 
the same conclusions on thermodynamic equilibrium apply. However, we 
note that, in this approach e.g., particles ratios involving particles with 
hidden strangeness are generally not well predicted, again implying non- 
equilibrium behavior. 

2.1. Statistical approach - grand canonical formalism 

The basic quantity required to compute the thermal composition of particle 
yields measured in heavy ion collisions is the partition function Z(T, V). In 
the Grand Canonical (GC) ensemble, 



where H is the Hamiltonian of the system, Qi are the conserved charges and 
fj,Q i are the chemical potentials that guarantee that the charges Qi are con- 
served on the average in the whole system. Finally (3 = 1 /T is the inverse 
temperature. The Hamiltonian is usually taken such as to describe a hadron 
resonance gas. For practical reasons, the hadron mass spectrum contains 
contributions from all mesons with masses below ~1.5 GeV and baryons 
with masses below ~2 GeV. In this mass range the hadronic spectrum is 
well established and the decay properties of resonances are reasonably well 
knowrPH This mass cut in the contribution of resonances to the partition 
function limits, however, the maximal temperature to r max ~ 200 MeV, up 
to which the model predictions may be considered trustworth y^ ' I w l 4z! l 3g l 




For higher temperatures the contributions of ( in general poorly known) 
heavier resonances are not negligible. The interaction of hadrons and reso- 
nances are usually only included by implementing a hard core repulsions, 
i.e. a Van der Waals-type interaction. Details of such a implementation are 
discussed below. The main motivation of using the Hamiltonian of a hadron 
resonance gas in the partition function is that it contains all relevant degrees 



Z GC (T, V, hq) = Tr[ e - p{H -^ 



(3) 
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off freedom of the confined, strongly interacting medium and implicitly in- 
cludes interactions that result in resonance formation. Secondly, this model 
is consistent with t he equ ation of state obtained from the LGT below the 
critical temperaturdZHEni J n a strongly interacting medium, one includes 
the conservation of electric charge, baryon number and strangeness. The 
GC partition function © of a hadron resonance gas can then be written 
as a sum of partition functions In Zi of all hadrons and resonances 

In Z(T, V,fl)=J2 ln Z ^ V > ( 4 ) 

i 

where e, = \Jp 2 + mf and fl = (/Ab, //s, //q) with the chemical potentials 
Hi related to baryon number, strangeness and electric charge, respectively. 
For particle i of strangeness Si, baryon number Bi, electric charge Qi and 
spin-isospin degeneracy factor 



ln Zi(T, V,p) = ^J ±p 2 dp ln[l ± A. ( exp(-/? ej 
with (+) for fernhons, (-) for bosons and fugacity 



(5) 



, frr BjHB + S^s + Q^q 

Xi(T, fi) = exp( 2.) (6) 

Expanding the logarithm and performing the momentum integration in 
Eq. (jHJl we obtain 

mz^v,® = ^E^a^ 2 (^), (7) 

fe=l 

where Ki is the modified Bessel function and the upper sign is for bosons 
and lower for fermions. The first term in Eq. (JJJ) corresponds to the Boltz- 
mann approximation. The density of particle i is obtained from Eq. Q 
as 

, ^ (m Tg. l ^(±l) k+1 . , ,km is 

ni (T,fi) = ^D^-W'fy). ( 8 ) 
fc=i 

The partition function @ is the basic quantity that allows to describe 
all thermodynamical properties of a fireball composed of hadrons and reso- 
nances being in thermal and chemical equilibrium. In view of further appli- 
cation of this statistical operator to the description of particle production 
in heavy ion collisions we write explicitly the results for particle density 
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Fig. 3. The ratio of the total density of positively charged pions that includes all 
resonance contributions to the density of thermal pions. The calculations are done in the 
hadron resonance gas model for fig =250, 550 MeV and for different temperatures. 

obtained from Eq. @ . Of particular importance here is to account for res- 
onances and their decay into lighter particles. The average number (Ni) 
of particles i in volume V and temperature T, that carries strangeness Si, 
baryon number Bi, and electric charge Qi, is obtained from Eq. as 

(AQ(T,£) = (N t ) th (T,t) + J2 Tj^iNjY^iT^) (9) 

where the first term describes the thermal average number of particles of 
species i and second term describes overall resonance contributions to par- 
ticle multiplicity of species i. This term is taken as a sum of all resonances 
that decay into particle i. The I^—^ is the corresponding decay branching 
ratio of j — ► i. The corresponding multiplicities in Eq. J!|J) are obtained from 
Eq. (jHJ). The importance of the resonance contribution to the total particle 
yield in Eq. © is illustrated in Fig. J2Jl as the ratio of total to thermal num- 
ber of 7r + . From this figure it is clear that at high temperature (or density) 
the overall multiplicity of light hadrons is indeed dominated by resonance 
decays. In the high-density regime, that is for large T and/or \ib, the repul- 
sive interactions of hadrons should be included in the partition function Q . 
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To incorporate the repulsion at short distances one usually uses a hard core 
description by implementing ex clud ed volume correctiond^l In a thermo- 
dynamically consistent approacfP2 these corrections lead to a shift of the 
baryon-chemical potential. . We discuss below how this is implemented in 
our calculations. The repulsive interactions are important when discussing 
observables of density type. Particle density ratios, however, are only weakly 
affectedSHI by the repulsive corrections. The partition function (@J depends 
in general on five parameters. However, only three are independent, since 
the isospin asymmetry in the initial state fixes the charge chemical poten- 
tial and the strangeness neutrality condition eliminates the strange chem- 
ical potential. Thus, on the level of particle multiplicity ratios we are only 
left with temperature T and baryon chemical potential as independent 
parameters. In Fig. J3J we show the relation of fj,s = ^s(T,hb) obtained 
from the strangeness neutrality condition. For low temperature this rela- 
tion is highly non-linear. For larger T, however, us shows an almost linear 
dependence on /is- One sees by inspection of Fig. (@J that, at T ~ 200 
MeV and hb ~ 300 MeV, /is ~ \^b- This relation is obtained in a QGP 
from strangeness neutrality conditions. In the present context of a hadron 
resonance gas this is a pure accident with no dynamical information. 

At lower energies, in pra ctise fo r T < 100 MeV, the widths of the 
resonances have to be includecP^ESl m Eq JgJ This is because the number 
of light particles coming from the decay of resonances is increased by the 
finite resonance width. In practice, the width of the A resonance is most 
importantP^EZI Thus, the approximation of the resonance width by a 5 
function is not justified. Assuming the validity of Boltzmann statistics one 
replaces the partition function in equation (JJJ by: 

In Z R = N — ^ T exp[(Sfl^ B + QrHq + S R fi s )/T] 

f dssK 2 (VS/T) - 7fj* (10) 

Jsmin ^ ( S ~ m R> + m R T R 

where s m i n is chosen to be the threshold value for the resonance decay and 
^/Smax ~ ™_r + 2Fr. The normalization constant iV is adjusted such that 
the integral over the Breit-Wigner factor gives 1. 
The statistical model, outlined above, 

was a ppiie oP-CHUBl 

to describe 

particle yields in heavy ion collisions. The model was compared with all 
available experimental data obtained in the energy range from AGS up to 
RHIC energy. Hadron multiplicities ranging from pions to omega baryons 
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7=770 MeV 




0.0 0.2 0.4 0.6 
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Fig. 4. The strange chemical potential us as a function of baryon— chemical potential 
for T=120,170 and 200 MeV. The results are obtained by imposing the strangeness 
neutrality condition in a hadron resonance gas. 

and their ratios were used to verify that there is a set of thermal param- 
eters (T,/j,b) which simultaneously reproduces all measured yields. In the 
following Section we present the most recent analysis of particle production 
in A-A collisions at RHIC, SPS and AGS energies. 



2.2. Thermal analysis of particle yields from AGS to 
RHIC energies 

For the analysis of data in the energy range of 40 GeV/nucleon and up- 
wards' 1 we use a grand canonical ensemble to describe the partition function 
and hence the density of the particles following Eqs. |Q]-|§J). As discussed 
above the temperature T and the baryochemical potential /is are the two 
independent parameters of the model, while the volume of the fireball V, 
the strangeness chemical potential fig, and the charge chemical potential /iq 
are fixed by the following additional conditions. First, overall strangeness 

d The results in this section were obtained in collaboration with D. Magestro and are 
published in part in Refs. 1351 1591 . 
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conservation fixes /xg. Note that this applies strictly for data integrated 
over 47T. For slices near mid-rapidity this condition is, however, also ap- 
propriate as the flow of strangeness in and out of the rapidity slice under 
consideration very nearly cancels. Charge conservation implies a condition 
on I 3 according to: 



Here, Z and N are the proton and neutron numbers of the colliding nuclei, 
I 3 and if are the third component of the total isospin and that of particle i. 
This condition is appropriate (and relevant) at lower beam energies where 
there is full stopping and An yields are used. For details see Refs. (I351lfif])l . At 
higher energies and for data analyzed in rapidity slices the right hand side 
of Eq. has to be replaced by the neutron excess of baryons transported 
into the rapidity slice under consideration. This number is clearly smaller 
than the full neutron excess entering Eq. but in general not well known. 
However, its precise knowledge is less relevant for higher beam energies since 
the isospin balance is dominated by pions. For practical purposes isospin 
conservation is important for AGS energies and below but its effect is small 
(on the 10 % level) already at 40 GeV/nucleon beam energy (where we have 
used as an upper limit the full neutron excess of the colliding nuclei, leading 
to a slight overestimate of the pion charge asymmetry) and negligible at top 
SPS and RHIC energies. Finally, the volume (which drops out anyway for 
particle ratios) can be obtained from total baryon number conservation (for 
full stopping and quantities which are evaluated over the complete phase 
space) or is fixed by using the measured pion multiplicity in the rapidity 
slice under consideration. As discussed above the hadronic mass spectrum 
used in the calculations extends over all mesons with masses below 1.5 
GeV and baryons with masses below 2 GeV. To take into account a more 
realistic equation of state we incorporate the repulsive interaction at short 
distances between hadrons by means of the excluded volume correction 
discussed above. A number of different corrections have been discussed in 
the literature. Here we choose that proposed in Refs. IjfiOl Ifill l%2"|) : 



p excL (T,n) = p ld - gas (T, fi); with H = f i-v eigen p excL (T,v). (12) 



This thermodynamically consistent approach to simulate interactions be- 
tween particles by assigning an eigenvolume v e i gen to all particles modifies 
the pressure p within the fireball. Equation (|12|l is recursive, as it uses the 




(11) 
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Fig. 5. Comparison between thermal model predictions and experimental particle 
ratios for Pb-Pb collisions at 40 GeV/nucleon. The thermal model calculations are ob- 
tained with T = 148 MeV and \x B = 400 MeV. 



modified chemical potential £i to calculate the pressure, while this pressure 
is also used in the modified chemical potential, and the final value is found 
by iteration. Particle densities are calculated by substituting n in Eq. (JSJ 
by the modified chemical potential /}. The eigenvolume has to be chosen ap- 
propriately to simulate the repulsive interactions between hadrons, and we 
have investigated the consequences for a wide range of parameters for this 
eigenvolume in Ref. 1)381 160f) . Note that the eigenvolume is v e i gen — 4^nR 3 
for a hadron with radius R. Assigning the same eigenvolume to all particles 
can reduce particle densities drastically but hardly influences particle ratios. 
Ratios may differ strongly, however, if different values for the eigenvolume 
are used for different particle species. Our approach here is, to determine, for 
nucleons, the eigenvolume acc ording to the hard-core volume known from 
nucleon-nucleon scattering^! Consequently, we assigned 0.3 fm as radius 
for all baryons. For mesons we expect the eigenvolume not to exceed that of 
baryons. For lack of better theoretical guidance we chose also for the mesons 
a radius of 0.3 fm. For a discussion of the implications of varying these ra- 
dius parameters see Ref. (fTTHl Ifffi)) . After thermal "production", resonances 
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Fig. 6. Comparison between thermal model predictions and experimental particle 
ratios for Pb-Pb collisions at 158 GeV/nucleon. The thermal model calculations are 
obtained with T = 170 MeV and fi B = 255 MeV. 



and heavier particles are allowed to decay, therefore contributing to the fi- 
nal particle yield of lighter mesons and baryons, as indicated above. Decay 
cascades, where particles decay in several steps, are also included. System- 
atic parameters regulate the amount of decay products resulting from weak 
decays. This allows to simulate the different reconstruction efficiencies for 
particles from weak decays in different experiments. In the following we 
compare predictions of the model with results of measured particle ratios 
for central Pb-Pb collisions at SPS energies (40 and 158 GeV/nucleon) and 
for central Au-Au collisions at RHIC energies Js nn = 130 and 200 GeV. 
An important issue in this context is whether to use data at mid-rapidity 
or data integrated over the full phase space. While it is clear that full 4ir 
yields should be used at low beam energies, this is not appropriate any more 
as soon as fragmentation and central regions can be distinguished. In that 
case the aim is to identify a boost-invariant region near mid-rapidity and to 
choose a slice in rapidity within that region. For RHIC energies this implies 
that an appropriate choice, given the available data, is a rapidity interval 
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of width Ay = 1 centered at midrapidity. The anti-proton/proton ratio 
stays essentially constant within that interval, b ut d rops rather strongly 
for larger rapidities and similar results are observedS2l for other ratios. Fur- 
thermore, t he r apidity distribution exhibits a boost-invariant plateau near 
mid- rapidityEl. As has been demonstrated in Ref. (65), effects of hydro- 
dynamic flow cancel out in particle ratios under such conditions. At SPS 
energies a boost-invariant plateau is not fully developed but stopping is not 
complete, either. In addition, the proton and anti-proton rapidity distribu- 
tions differ rather drastically, especially near the fragmentation regions, 
implying that particle ratios depend on rapidity (see. e.g., Ref. (66)). Un- 
der those circumstances we have decided to use, wherever available, data 
in a slice of ±1 unit of rapidity centered at mid-rapidity. This is slightly 
different from the analysis performed in Ref. I)38|l . where both mid-rapidity 
and fully integrated data were used. We note, however, see below, that the 
fit parameters T and /is obtained at 158 A GeV are very close to those 
determined earlier. The criterion for the best fit of the model to data was 
a minimum in 



culated from our model or measured in the experiment, and Oj represent 
the errors (including systematic errors where available) in the experimen- 
tal data points as quoted in the experimental publications. For the data 
we used all information available including that presented at the QM2002 
conference in July 2002. Details on the data selection, corrections for the 
weak-decay reconstruction efficiency, as well relevant references are found 
in Ref. (|59[) . Under the conditions discussed above the data can all be well 
described, as is detailed below, by a thermal distribution with T and [Ib 
as independent parameters. There is no need to introduce additional pa- 
rameters such as as strangeness suppression factors. The results of the fits 
for central Pb-Pb collisions at 40 and 158 GeV per nucleon are presented 
in Figs. I|5I6|) . At 40 GeV/nucleon 11 particle ratios are included in the fit, 
while the number is 24 at 158 GeV/nucleon. We obtain values for (T,/j,b) 
of (148±5, 400±10) and (170±5, 255±10), respectively, with reduced \ 2 
values of 1.1 and 2.0. Obviously the fits are quite good. A possible excep- 
tion is the <t>/ (ir + +ir~) ratio at top SPS energy, where there are conflicting 




(13) 



In the above equations 7^.™ odcl and IZ^ 



are the ith particle ratio as cal- 
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Fig. 7. Comparison between thermal model predictionJSH and experimental particle 
ratios for Pb-Pb collisions at ^/s nn = 130 GeV. Calculations were performed for T = 174 
MeV and fi B = 46 MeV. 



data from NA49 and NA50. This is already discussed in detail in Ref. I)38|l 
and no new information on this problem has appeared since. Note that this 
ratio has not been used in the x 2 minimization. The somewhat larger val- 
ues of x 2 at full SPS energy and the remaining uncertainty in fxg are due 
to a systematic problem not yet sufficiently addressed by the experiments. 
The contribution of weak decays of strange baryons to final baryons has 
been discussed and cuts are applied to reduce this contribution in the data. 
However, there is also a contribution from weak decays to charged pions 
(or, more generally) charged hadrons which is up to now poorly quantified 
by the experiments. If, e.g., in the ratio A/h~ , feeding of the A by decays of 
El's is suppressed due to cuts, but the n~ measurement has a 50 % efficiency 
for detection of pions from weak decays, the ratio would drop by 15-20 %, 
compared to the case with % eficiency for weak decays. With this option 
the reduced \ 2 value for the 158 GeV fit would drop from 2.0 to 1.5. This 
discussion indicates that there are sources of systematic uncertainties not 
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included in the data. The corrections for weak decays are, consequently, of 
utmost importance when discussing the "precision" of fits. 
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Fig. 8. Comparison of the experimental data on different particle multiplicity ratios 
obtained at RHIC at \/*jvjv = 130 and 200 GeV with thermal model calculations. The 
thermal model analysis is from Refs. J3S]l2fif an d recent update by D. Magcstro. 

The results for RHIC energies are shown in Figs. |Q|SJ|. In Fig. Q we 
present the results as published in Ref. I|35[) in the summer of 2001. Since 
then, the data at ^/s nn = 130 GeV have been consolidated and extended 
and first (in some cases still preliminary) results have been provided for 
y^Snn = 200 GeV. The current state of affairs in summarized in Fig. J7J • The 
results demonstrate quantitatively the high degree of equilibration achieved 
for hadron production in central Au-Au collisions at RHIC energies. We 
obtain values for (T> B ) of (174±7, 46±5) and (177±7, 29±6), respectively, 
with reduced \ 2 values of 0.8 and 1.1. We note that ratios involving multi- 
strange baryons are well reproduced as is the <fi/h~ ratio. Even relatively 
wide resonances such as the K*'s fit well into the picture of chemical freeze- 
out. This obviates the need for quark coalescence models as proposed in 
Ref. ljr)?j> and non-equilibrium models as proposed in Ref. (68). 

Very recently, the STAR collaboration has provider^ first data, with 
about 30 - 50 % accuracy, on the p°/tt and /°(980)/tt ratios in semi-central 
Au-Au collisions. These mesons have been reconstructed in STAR via their 
decay channel in 2 charged pions. Comparing the preliminary results from 
STAR with our thermal model prediction reveals that the measured ratios 
exceed the calculated values by about a factor of 2. This is quite surprising, 
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especially considering that we use a chemical freeze-out temperature of 177 
MeV for the calculation, while one might expect these wide resonances 
to be formed near to thermal freeze-out, i.e. at a temperature of about 
120 MeV. At this temperature, the equilibrium value for the p°/tt ratio is 
about 4 • 1CP 4 , while it is 0.11 at 177 MeV. Even with a chemical potential 
for pions of close to the pion mass and taking into account the apparent 
(downwards) mass shift of 60 - 70 MeV for the p° it seems difficult to 
explain the experimentally observed value of about 0.2. 

We finally note that the model discussed here was also applied to the 
AGS data collected in Ref. I|37(l . The best fit, obtained for Rbaryon=Rmeson- 
=0.3 fm, yields T = 125 (+3-6) MeV and [i B = 540 ±7 MeV, well in line 
with the calculations reported in Ref. (37.) ■ 

In summary, hadron multiplicities produced in central nucleus-nucleus 
collisions in the range of AGS to full RHIC energy can be quantitatively 
described with a grand-canonical partition function based on the full hadron 
resonance spectrum, assuming complete chemical equilibrium. There is no 
need to introduce non-equilibrium parameters or strangeness suppression 
factors if data near mid-rapidity are considered. The physical relevance of 
the two model parameters T and /is is described in detail in our discussions 
below concerning the phase boundary between hadrons and the quark-gluon 
plasma. 

2.3. Comparison of measured particle densities with 
thermal model predictions 

As discussed below, the value for the energy density predicted by the 
presently used thermal model, including the excluded volume correction, 
agrees well with results from the lattice for temperatures below the critical 
temperature. It makes therefore sense to compare the densities for pions and 
nucleons predicted by the model with values determined from experiments. 
The CERES collaboration has recently performed an analysis of 2-pion cor- 
relation experiments for the energy range between AGS a nd RHIC, from 
which values for these densities have been determine J 7U I 71 I from data taken 
at mid-rapidity. For the nucleon density (at thermal freeze-out) the exper- 
imental numbers are, at 40 and 158 GeV/nucleon c and at ^Js nn = 130 



c Wc take here the data published in Ref. 1661 : the data reported in Ref. 1721 are about 
20 % lower and would not fit the beam energy systcmatics. 
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GeV, 0.077 ±0.005/fm 3 , 0.063 ±0.005/fm 3 and 0.06 ±0.009/fm 3 . From the 
model we deduce, at chemical freeze-out, values of 0.10/fm 3 , 0.10/fm 3 and 
0.08/fm 3 . This would imply a volume increase of about 40 % from chemical 
to thermal freeze-out. For pions the situation could be more complicated 
since yield ratios involving pions are (apparently) fixed at chemical freeze- 
out, implying the build-up of a pion chemical potential between chemical 
and thermal freeze-out. From the data one 

deduce!!^ 

a pion density at 

(thermal) freeze-out of 0.28 ±0.03/fm 3 , 0.43 ±0.03/fm 3 , and 0.49 ±0.1/fm 3 , 
at 40 and 158 GeV/nucleon and at y/s nn = 130 GeV. These values should 
be contrasted with the calculated (chemical) freeze-out values of 0.35/fm 3 , 
0.59/fm 3 and 0.62/fm 3 . From these numbers one would conclude a 30 % 
volume increase between chemical and thermal freeze-out, assuming that 
the pion chemical potential fixes the pion number to the value obtained 
at chemical freeze-out. This rather small volume increase indicates that 
the time between chemical and thermal freeze-out cannot be very long at 
SPS and RHIC energies. At AGS energy, the corresponding 7r + and proton 
densities of 0.051/fm 3 and 0.053/fm 3 agree well with those estimatec F^^ 
from particle interferometry (0.058/fm 3 and 0.063/fm 3 , respectively) im- 
plying that, at AGS energy, thermal and chemical freeze-out take place at 
nearly identical times and temperatures. 



2.4. Statistical model and composite particles 

An often overlooked aspect of the thermal model is the possibility to com- 
pute also the yields of composite particles. For example, the d/p and d/p 
ratios measured at SPS and AGS energies are we ll reproduced^ with the 
same parameters which are used to describd^^ baryon and meson ratios. 
Furthermore, the AGS E864 Collaboration has recently published*^ yields 
for composite particles (light nuclei up to mass number 7) produced in cen- 
tral Au-Au collisions at AGS energy near mid-rapidity and at small pt ■ In 
this investigation, an exponential decrease of composite particle yield with 
mass is observed over 7-8 order of magnitude, yielding a penalty factor P p 
of about 48 for each additional nucleon. Extrapolation of the data to large 
transverse momentum values, considering the observed mass dependence 
of the slope constants, reduces this penalty factor to about 26, principally 
because of transverse flow. In the thermal model, this penalty factor can 
be related with thermal particle phase-space. In the relevant Boltzmann 
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approximation, we obtain 

m±fi b 

-Rp-exp — - — , (14) 

where m is the nucleon mass and the negative sign applies for matter, 
the positive for anti-matter. Small corrections due to the spin degeneracy 
and the A 3 / 2 term in front of the exponential in the Boltzmann formula for 
particle density are neglected. Using the freeze-out parameters T=125 MeV 
and fib = 540 MeV appropriate^ for AGS energy one getJ^3 Rp ~ 23, in 
close agreement with the data for the produ ctio n of light nuclei. It was also 
noted that the anti-matter yields measurecpSI by the E864 Collaboration 
yield penalty factors of about 2-10 5 , again close to the predicted^ value of 
1.3-10 5 . 

This rather satisfactory quantitative agreement between measured rel- 
ative yields for composite particles and thermal model predictions provides 
some confidence in the predictions for yields of exotic objects produced in 
central nuclear collisions. We briefly comment here on the results obtained 
in Ref. <E3). 

In this investigation, the production probabilities for exotic strange ob- 
jects and, in particular, for strangelets were computed in the thermal model. 
The results are reproduced in Table ^ for temperatures relevant for beam 
energies between 10 and 40 GeV/nucleon. We first note that predictions of 
the thermal model and, where available, the coalescence model of Ref. I|98|l 
agree (maybe surprisingly) well particularly for lighter clusters. Secondly, 
inspection of Tabled also shows that, in future high statistics experiments 
which will be possible at the plannec™ new GSI facility, multi-strange 
objects such as ^AA^e should be experimentally accessible with a planned 
sensitivity of about 10~ 13 per central collision in a years running, should 
they exist and be produced with thermal yields. Investigation of yields of 
even the lightest conceivable strangelets will be difficult, though. 



3. Exact Implementation of the conservation laws in the 
statistical models 

The analysis of particle yields obtained in central heavy ion collisions from 
AGS up to LHC energy has shown that hadron multiplicities are very well 
described by assuming a complete thermalized state at fixed T and /is ■ In 
this broad energy range, particle yields and their ratios are, within experi- 
mental error, well reproduced by the statistical hadron resonance gas model 
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Table 1. Produced number of nonstrange and strange clusters and of strange quark 
matter per central Au+Au collision at AGS energy, calculated in a thermal model for 
two different temperatures, baryon chemical potential fii> = 0.54 GeV and strangeness 
chemical potential fi s such that overall strangeness is conserved. 

Thermal Model Parameters 



Particles 
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Source: The Coalescence model predictions in the last column are from Table 2 of 
Ref. (HHJ- 

that accounts for the conservation laws of baryon number, strangeness and 
electric charge in the grand canonical ensemble. The natural question aris- 
ing here is whether this statistical order is a unique feature of high energy 
central heavy ion collisions or is it also there at lower energies as well as in 
hadron-hadron and peripheral heavy ion collisions. To address this question 
one needs, however, to stress that when going beyond high energy central 
heavy ion collisions the grand canonical statistical operator 10 has to be 
modified. 

Within the statistical approach, particle production can only be de- 
scribed using the grand canonical ensemble with respect to conservation 
laws, if the number of produced particles that carry a conserved charge is 
sufficiently large. In view of the experimental data this also means that 
the event-averaged multiplicities are controlled by the chemical potentials. 
In this description the net value of a given charge (e.g. electric charge, 
baryon number, strangeness, charm, etc.) fluctuates from event to event. 
These fluctuations can be neglected (relative to the squared mean parti- 
cle multiplicity) only if the particles carrying the charges in question are 



24 



Peter Braun-Munzinger, Krzysztof Redlich, Johanna Stachel 



abundant. Here, the charge will be conserved on the average and the grand 
canonical description developed in the last section is adequate. In the op- 
posite limit of low production yield the particle number fluctuation can be 
as large as its event averaged value. I n this c ase charge conservation has to 
be implemented exactly in each even 

The exact conservation of quantum numbers introduces a constraint on 
the thermodynamical system. Consequently, the time dependence and equi- 
librium distribution of particle multiplicity can differ from that expected 
in the grand canonical limit. To see these differences one needs to per- 
form a detailed study of particle equilibration in a thermal environment. 
To discuss equilibration from the theoretical point of view one needs to 
formulate the kinetic equations for particle production and evolution. In a 

partonic medi um th is requires, in general, the formulation of a transport 
equation" »^<"H "4 

involving colour degrees of freedom and a non-Abelian 
structure of QCD dynamics . In the hadronic medium, on the other hand, 
one need J^* ^ W l H?l | H?5 1 1 071 to account for the charge conservations related 
with the U(I) internal symmetry. 



3.1. Kinetics of time evolution and equilibration of charged 
particles 

In this section we will discuss and formulate the kinetic equations that 
include constraints imposed by the conservation laws of Abelian charges 
related with U(l) internal symmetry. We will indicate the importance of 
the conservation laws for the time evolution and chemical equilibration 
of produced particles and their probability distributions. In particular, we 
demonstrate that the constraints imposed by the charge conservation are of 
crucial importance for rarely produced particle species such as for particles 
with hidden quantum numbers like e.g. for J/ip. 

To study chemical equilibration in a hadronic medium we introduce 
first a kinetic model that takes into account the production and annihila- 
tion of particle-antiparticle pairs cc carrying U(l) quantum numbers like 
strangeness or charm. It is also assumed that particles c and c are produced 
according to a binary process ab — > cc and that all particle momentum 
distributions are thermal and described by the Boltzmann statistics. The 
charge neutral particles a and b are constituents of a thermal fireball with 
temperature T and volume V. We will consider the time evolution and 
equilibration of particles c and c inside this fireball, taking into account the 
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constraints imposed by the U(l) symmetry. First, we formulate a general 
master equation for the probability distribution of particle multiplicity in 
a medium with vanishing net charge and consider its properties and solu- 
tions. Then we will discuss two limiting cases of abundant and rare particle 
production. Finally, the rate equation will be extended to a more interesting 
situation where there are different particle species carrying the conserved 
quantum numbers inside a thermal fireball that also has a non vanishing 
net charge. 

3.1.1. Kinetic master equation for probabilities 

Consider P/v c (r) as the probability to find N c particles c, where < N c < 
oo. This probability will obviously change in time owing to the production 
ab — > cc and absorption cc — > ab processes. The equation for the probability 
Pjv c contains terms which increase in time, following the transition from 
N c — 1 and N c + 1 states to the N c state, as well as terms which decrease 
since the state N c can make transitions to N c + 1 and N c — 1 (see Fig.EJ). 




Fig. 9. A schematic view of the master equation for the probability Pjv(t) due to 
ab <-> cc and the inverse process. 

The rate equation is determined by the magnitude of the transition 
probability per unit time due to the production G/V and the absorption 
L/V of cc pairs through ab <-> cc process. The gain (G =< o- Q (,^ccW a fc >) 
and the loss {L =< o~ C c-^abVcc >) terms represent the momentum average 
of particle production and absorption cross sections. 

The transition probability per unit time from N c + 1 — > N c is given by 
the product of the probability L/V that the single reaction cc — > ab takes 
place multiplied by the number of possible reactions which is formally, 
(N c + 1)(N S + 1). In the case when the charge carried by particles c and 
c is exactly and locally conserved, that is if ( N c + N s = 0) , this factor is 
just (N c + l) 2 . Similarly, the transition probability from N c — > N c + 1 is 
described by G(N a )(Nb)/V, where one assumes that particles a and b are 
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not correlated and their multiplicity is governed by the thermal averages. 
One also assumes that the multiplicity of a and b is not affected by the ab — > 
cc process. The mast er e quation for the time evolution of the probability 
Pjv c (r) can be writterPSin the following form: 

^ = ^(N a )(N b )P N ^ + ±{N C + l) 2 P Ne+1 

- y{N a )(N b )P Nc ~^N*P Nc . (15) 

The first two terms in Eq. I|15|) describe the increase of Pn c (t) due to the 
transition from N c — 1 and N c + 1 to the N c state. The last two terms, on 
the other hand, represent the decrease of the probability function due to 
the transition from iV c to the N c + 1 and N c — 1 states, respectively. 

For a thermal particle momentum distribution and under the Boltzm ann 
approximation the thermal averaged cross sections are obrainecP^^ from 



< <r ab ^ccv ab >- 8 mlmlK 2 ((3m a )K 2 (f3m b ) ' (ibj 

where K\, K 2 are modified Bessel functions of the second kind, m^ b = 
m a + m b and m~ b = m a — m b , t = s/s is the center-of-mass energy, 
P the inverse temperature, v ab = {{k a k b ) 2 — m^ml)/E a E b is the rela- 
tive velocity of incoming particles and the integration limit is taken to 
be to = max[{m a + m b ), (m c + m e )]. 

The rate equation for probabilities l|15|) provides the basis to calculate 
the time evolution of the momentum averages of particle multiplicities and 
their arbitrary moments. Indeed, multiplying the above equation by iV c 
and summing over N c , one obtains the general kinetic equation for the 
time evolution of the average number (N c ) — J^n =o N c Pn c (t) of particles 
c in a system. This equation reads: 

^ L = y(N a )(N b )-^ c ). (17) 

The above equation cannot be solved analytically as it connects particle 
multiplicity (N c ) with its second moment (N%). However, solutions can be 
obtained in two limiting situations: i) for an abundant production of c 
particles, that is when (N c ) ^ 1 or ii) in the opposite limit of rare particle 
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production corresponding to (N c ) <C 1. Indeed, since 

{N*) = {N C ) 2 + {6N*), (18) 

where (SN 2 ) represents the fluctuations of the number of particles c, one 
can make the following approximations: 

i) for (N c ) > 1 one has a (JV|) w (iV c ) 2 , and Eq. JEJ obviously reduces 
to the well knowif^form: 

(19) 

ii) however, for the rare production, particles c and c are strongly correlated 
and thus, for (N c ) <C 1 one takes (N%) ~ (-Nc)j consequently Eq. I|17|l takes 
the form: 

d(N c ) G , . . . L . 

^«-(7V a )(7V b )--(iV c ), (20) 

where the absorption term depends only linearly, instead of quadratically, 
on the particle multiplicity. 

From the above it is thus clear that, depending on the thermal conditions 
in the system (that is its volume and temperature) , we are getting different 
results for the equilibrium solution and the time evolution of the number 
of produced particles c. This is very transparent when solving the rate 
equations lfT3jl and J2DJ. 

In the limit wh en (iV r ) 3> 1, the standard Eq. Ijl9|l is valid and has the 
well known solutio: 



(N c (t)) = (7V c ) cq tanh (r/r ), (21) 

where the equilibrium value (N c ) cq of the number of particles c and the 
relaxation time constant r are given by: 

(N c ) cq = V~e , T ° = Yj= e > ( 22 ) 

respectively, with e = G(N a )(Nb)/L. 

In the particular case when the particle momentum distribution is ther- 
mal, the ratio of the gain (G) to the loss (L) terms can be obtained from 
Eq. JH| as 

G = d c m 2 c K 2 (m c /T)d E m 2 5 K 2 (m s /T) 
L d a mlK 2 {m a /T)d b mlK 2 (m b /Ty 
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where we have employed the detailed balance relation between the cross 
sections for production a a b and for absorbtion a cS for ab «-> cc processes 

( v _ dgrfb [t 2 - (m+ 5 ) 2 ]^ 2 - (?71~ e ) 2 ] 
dcrfe [i 2 - (m+ ) 2 ][t 2 - (m afc ) 2 ] 

with <ii being the spin-isospin degeneracy factor and as in Eq. (|16|l . 

In Boltzmann approximation, the equilibrium average number of parti- 
cles c in Eq. 122(1 reads: 

Weq = ^T m ^ 2 (m c /T). (25) 

This is a well known result for the average number of particles in the Grand 
Canonical (GC) ensemble with respect to the U(l) internal symmetry of 
the Hamiltonian. The chemical potential, which is usually present in the 
GC ensemble, vanishes in this case, because of the requirement of charge 
neutrality of the system. Thus, the solution of Eq. fB^I results in the ex- 
pected value for the equilibrium limit in the GC formalism where a charge 
is conserved on the average. 

In the opposite limit, where (N c ) <C 1, the time evolution of a particle 
abundance is described by Eq. (121 that has the following solution: 

(7V c (r)) c = (iV c ) c c q (l-e-^), (26) 

with the equilibrium value and relaxation time given by 

<Ag c C q = e, T o=J- (27) 

The above result, as will be shown in the next section, is the asymptotic 
limit of the particle multiplicity obtained in the canonical (C) formulation 
of the conservation lawj™. Here the charge related with the U(l) sym- 
metry is exactly and locally conserved, contrary to the GC formulation 
where this conservation is only valid on the average. 

Comparing Eq. with Eq. we first find that, for (iV c ) < 1, 

the equilibrium value is by far smaller than what is obtained in the grand 
canonical limit, i.e. 



(iV c ) e c q = (N c ) cq 2 « (N c ) cq . (28) 
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Secondly, we can conclude that the relaxation time for a canonical sys- 
tem is shorter than the grand canonical value, i.e. 

r ° = T-oWc q « t , (29) 

since in the limit (ii) the equilibrium value (N c ) eq <C 1. 

We note that the (i) and the (ii) limits are essentially determined by 
the size of (SN%), the fluctuations of the number of particles c. The grand 
canonical results correspond to small fluctuations, i.e. (8N 2 ) / (N c ) 2 < I, 
while large fluctuations (8N 2 ) / (N c ) 2 > 1 require a canonical description. 

The volume dependence of particle density obviously differs in the C 
and in the GC limit. The particle density in the GC limit is ^/-independent 
whereas in the canonical approach it can even scale linearly with V. 

The difference between the C and asymptotic GC result already seen 
on the level of the rate equations (|17I19() . is even more transparent when 
comparing master equations for probabilities. In the following we formulate 
this equation for the GC description of quantum number conservation. 

In case of abundantly produced particles c and c through the ab — > cc 
process we do not need to worry about strong particle correlations due 
to charge conservation. This also means that, instead of imposing charge 
neutrality conditions through N c — N E = 0, one assumes conservation on 
the average, that is (N c ) — (N s ) = 0. In this case the master equation (fT5|) 
can be simplified. 

In the derivation of Eq. Ijlfijl the absorption terms proportional to L were 
obtained by constraining the charge conservation to be local and exact. 
For the conservation on the average, the transition probability from N c 
to the (N c — I) state is no longer proportional to (L/V)N 2 but rather to 
(L/V)N C (N S ), since the exact conservation condition N c = N s is no longer 
valid and the number of c particles can only be determined by its average 
value. In the GC limit, the master equation for the time evolution of the 
probability P/v c (r) takes the following form: 



^ = ^(N a )(N b )P Nc ^ + ^(N c + l)(N 5 )P Nc+1 

-y(N a )(N b )P Nr ~l±N c (N 5 )P Nc . (30) 

Multiplying the above equation by N c , summing over N c and using the 
condition that (N c ) = {Nc}, one recovers Eq. lfT9")l. the rate equation for 
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(N c ) in the GC ensemble. The above equation is thus indeed the general 
master equation for the probability function in the GC limit. Comparing 
this equation with the more general Eq. (|15|l . one can see that the main 
difference is contained in the absorption terms that are linear in particle 
number instead of being quadratic . 

Eq. (|3U[1 can be solved exactly. Indeed, introducing the generating func- 
tion g(x,r) for Pn c , 



g(x,r)=J2 xNcp N a (r), (31) 

N c =0 

the iterative equation H3()|l for the probability can be converted into a dif- 
ferential equation for the generating function: 

^^ = ^(l-*)[</-^], (32) 
with the general solution^? 

9(x, t) = ffo (1 - xe~ T ) exp[\/e(l - x)(e~ T - 1)], (33) 

where g' = dg/dx, f = {Ly/e/V)r and ^fi = (N c ) eq given by Eq. (9). 

One can readily find out an equilibrium solution to the above equation. 
Taking the limit t = oo in the Eq. (J33J) leads to 

3 oq (a;) = exp[-Ve(l - x)}, (34) 

with the corresponding equilibrium multiplicity distribution: 



This is the expected Poisson distribution with average multiplicity ^fl 



PN a , eq = ^r^- (35) 



3.1.2. The equilibrium solution of the general rate equation 

The master equation 13U|) . that describes the evolution of the probability 
function in the GC limit, could be solved analytically. The general equation 
(|15|l . however, because of the quadratic dependence of the absorption terms, 
requires a numerical solution. Nevertheless, the equilibrium result for the 
particle multiplicity can be given. 
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Converting Eq. I|15|) for Pjf c into a partial differential equation for the 
generating function 



9 (x,t)=J2x n °P Nc (t). (36) 



N c =0 



one find; 



J29 



dg(x,r) L „ , 

— 7^ — = yQ--z)(xg +9 -m- (37) 
The equilibrium solution g eq (x) thus obeys the following equation: 

%9eq + .9cq - £ .9cq = 0. (38) 

By a substitution of variables (a; = y 2 e/4), this equation is reduced to the 
Bessel equation, with the following solution: 

9eq(x) = ^- Io{2y/ex), (39) 

where the normalization is fixed by g(l ) = ^ Pn c = 1< 

The equilibrium value for the probability function is now written 
from Eqs. as: 

= WWF" (40) 

We note that the equilibrium distribution of the particle multiplicity 
is not Poissonian. This fact was indicated first in equilibrium studies in 
Ref. <P,09). In our case this is a direct consequence of the quadratic depen- 
dence on the multiplicity in the loss terms of the master equation l|15|) . The 
Poisson distribution is obtained from Eq. (|40|l if ^/e 3> 1, that is for large 
particle multiplicity where the C ensemble coincides with the GC asymp- 
totic approximation. In Fig. i|l(J|) we compare the Poisson distribution from 
Eq. H35|) with the distribution from Eq. (|40|) for two values of y/i. 

The result for the equilibrium average number of particles c can be 
obtained as: 

w^-'v-S'Wffi ■ <41 » 

The above expression will be shown in the next section to coincide with 
the one expected for the particle multiplicity in the canonical ensemble with 
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Fig. 10. The probability function from Eqs. I35I4UI for two values of e = 4 and 16. 
The full lines represent Poisson distribution. 



respect to U(l) charge conservation I I The rate equation formulated 
in Eq. (|15f) is valid for arbitrary values of (N c ) and obviously reproduces 
(see Eqs. H101l - fTUl|l l the standard grand canonical result for a large (N c ). 
Thus, within the approach developed above one can study the chemical 
equilibration of charged particles following Eq. (JTSJ, independent of thermal 
conditions inside the system. 

3.1.3. The master equation in the presence of the net charge. 

So far, in constructing the evolution equation for probabilities, we have 
assumed that there is no net charge in the system under consideration. For 
the application of the statistical approach to particle production in heavy 
ion and hadron-hadron collisions, the above assumption has to be extended 
to the more general case of non- vanishing initial values of conserved charges. 
In the following we construct the evolution equation for Pjy (t) in a thermal 
medium assuming that its net charge S is non-vanishing. 

The presence of a non-zero net charge requires modification of the ab- 
sorption terms in Eq. I)15[l. The transition probability per unit time from 
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the N c to the N c — 1 state was proportional to (L/V)N C N S . Admitting 
an overall net charge 5^0 the exact charge conservation implies that 
N c — N s = S. The transition probability from N c to N c — 1 due to pair 
annihilation is thus (L/V)N C (N C — S). Following the same procedure as in 
Eq. (|15fl one can formulate the following master equation for the probability 
P^ (t) to find N c particles c in a thermal medium with a net charge S: 



dP s C J 
j± = -(N^iN^Pl^ + -(N c + 1)(N C + 1 - S)Pl +1 



^(N a )(N b )Pl -- Nc (N c -S)Pl, (42) 

which obviously reduces to Eq. I|15|) for S = 0. 

To get the equilibrium solution for the probability and multiplicity, we 
again convert the above equation to the differential form for the generating 
function g s (x, r) = J2n c =o xNcP n c ( t ) : 



= £(1 - x) (xg'£ + g' s (l -S)- eg s ) . (43) 

In equilibrium, dg (x,t)8t = and the solution for g^ q can be found as 
follows: 



S /2 

9c q (x) = j^ T) Is(2V^), (44) 

where the normalization is fixed by g(l) = Y] P n = 1. 

The master equation for the probability to find N 5 antiparticles c, its 
corresponding differential form and the equilibrium solution for the gener- 
ating function can be obtained by replacing S with — S in Eqs. (|42H44|) 

The result for the equilibrium average number of particles (N c ) eq and 
antiparticles (N s ) cq is obtained from the generating function using the re- 
lation: (N c ) cq — g'(l). The final expressions read: 



The charge conservation is explicitly seen by taking the difference of these 
equations that results in the net value of the charge S. 
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The thermal average values of the particle number given through Eq. 14511 
will be later deri ved from the equilibrium partition function by using the 
projection methocPSalllZI 



3.1.4. The kinetic equation for different particle species 

The rate equations discussed until now were derived assuming that there is 
only one kind of particle c and its antiparticle c that carry conserved charge. 
To study equilibration of particles in a strongly interacting environment one 
also needs to include processes that involve different species. In low energy 
heavy ion collisions e.g. the K + mesons are not only produced in pairs 
together with K~ but also with the strange hyperon A or E° due to the 
ttN — > AK + process. The contribution of K° and K° has to be also included 
as these particles are produced with similar strength as charged kaons. To 
account for this situation one generalizes the rate equations described in 
the last sections. 

Consider Pn k ,n a ( t ) as the probability to find Nk and N A number 
of K~ mesons and A baryons. Including the production and absorption 
processes such as: mm — > K + K~ and mN — > K + A this probability will 
obviously change in time. Here m(N) denotes a meson (nucleon). 
Following a similar procedure as was explained in Fig. © the master equa- 
tion for the time evolution of the probability Pm„ jv» (t) can be written 
assent)! 



= 9^ {Nm){Nfh)P + ^ {Nr + 1){Nk + l + N A ) 

UT V V 

- ^(N m )(N rll )P NK , NA - ^N K (N K + N a )P Nk>Na 

+ ^-{N m )(N N )P NxirfA -i + ^(N A + l)(N A + l + N K ) 
xPat k ,jv a +i (46) 

- ty-(N m )(N N )P NKtffA - ^-N A (N A + N k )P Nk , Na . 

with G m and L m being the production and absorption terms for the mm ^ 
K + K~ reaction and Gn and Lm denote equivalent terms for the mN # 
K + A process. 
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The equilibrium solution for the probability function Pn k n a can be 

Z? +Na (N K + N A )\e%<efr 



found as follows-^ 



e^ +NA (N K + N A y.£«e»- 

Nk ' Na /o(2^)((^+iV A )!) 2 e"f +NA N K \N A \ 



(47) 



with e to t — e m + £n and e m (N) — G m {N)(N mi ( m ))(N m2 ^ N - ) ) /L m ^ N y 

The equilibrium probability distribution is thus, according to Eq. 147fl . 

the product of the distribution of the number of pairs (Nk + N A ) and a 

binomial distribution that determines the relative weight of the individual 

particles, in our case the K~ and A. 

The probability Pjj is obviously normalized such that J^i j Pi.j = 1- 

The equilibrium value for the multiplicity (Ni) with i = K~ or i = A can 

be obtained as: 

\^K-)eq — —7= — T ln , n > (-^A/eq - — p — r , , x (48) 

with e m , c^v and defined as above. 

The average value of K + can be obtained applying strangeness conser- 
vation leading to: 



(N K+ ) eq = {N K -) eq + (N A ) eq . (49) 

The results presented here can be extendecP^l to an even more general 
case where there is an arbitrary number of different particle species carrying 
the quantum numbers related with U(l) symmetry of the Hamiltonian. 



3.2. The canonical description of an internal symmetry - 
projection method 

Using the above kinetic analysis of charged particle production probabil- 
ities we have demonstrated that equilibrium distributions does not neces- 
sarily coincide with the GC value. It is thus natural to ask what is the 
corresponding partition function that can reproduce the kinetic results ob- 
tained in Eqs. H41I45I48(1 . The main step in deriving these equations was 
an assumption of an exact conservation of quantum numbers in the kinetic 
master equations (|15I42I47|) . Thus, one should account for this important 
constraint in constructing the partition function. 
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The exact treat ment of q uantum numbers in statistical mechanic s has 
been well establishecPS^EEIfor some time now. It is in general obtainecP^ZEEl 
by projecting the partition function onto the desired values of the conserved 
charge by using group theoretical methods. In this section we develop these 
methods and show how one gets the partition function that accounts for 
exact conservation of quantum numbers. The derivation will be not only 
restricted to the charge conservation related with an Abelian U(l) internal 
symmetries and their direct products, but it will include also symmetries 
that are imposed by any semi-simple compact Lie group. 

The usual way of treating the problem of quantum number conserva- 
tion in statistical physics is by introducing the grand canonical partition 
function, as in Eq. (0. For only one conserved charge, e.g. strangeness S, 

Z(/! S ,T) = Tr[e-' 3( ^-^ ) ] (50) 

The chemical potential fis is then fixed by the condition that the average 
value of strangeness of a thermodynamical system is conserved and has the 
required value (S) such that: 

= T dlnZ(^T) 
0/J.S 

This method, as shown in the previous sections, is only adequate if the 
number of particles carrying strangeness is very large and their fluctuations 
can be neglected. 

In order to derive a partition function that is free from the above re- 
quirements let us first reorganize Eq. I|50(l . Denoting the states under the 
trace as \s) such that H\s) = E s \s) and Sjs) = s\s) one writes 

s— +oo S— +OC 

Z{n s ,T)= J2 e~^e^ = £ Z s \ s s (52) 

S— — QO S — — 00 

where we have introduced the fugacity A5 = e^^ s and where 

Z s = Tr s [e-P 6 ] (53) 

is just the partition function that is restricted to a specific total value S of 
the conserved charge. This is the canonical partition function with respect 
to strangeness conservation. Thus, Z$ is a coefficient in the Laurent series in 
the fugacity. Our goal is to calculate Zs- This is an easy task: starting from 
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Eq. (|52|l we apply the Cauchy formula and take an inverse transformation 
to obtain 

Zs(T ' v) = hj l^ z(As ' Tj v) (54) 

Choosing the integration path as the unit circle and parameterizing it as 
As = exp we can convert the contour integral into the angular one as 

Zs(T,V)= [ " p.Ztf,T,V) (55) 

J-7T ^ 

where the generating function Z((f>,T,V) = Z(\$ = e ic ^,T 7 V) is obtained 
from the grand canonical partition function by a Wick rotation of the chem- 
ical potential /is — > i<f>. This generating function is the same for all canon- 
ical partition functions with an arbitrary but fixed value of the conserved 
charge. Eq. Q55JI is the projection formula onto the canonical partition func- 
tion that accounts for the exact conservation of an Abelian charge. This is 
a projection procedure as Zg is obtained from 

Z S (T, V) = Tr s [e-P A ] = Tr [e~^ P s ] (56) 

where Ps = Pg is the projection operator on the states with the exact value 
of S. For an Abelian symmetry, Pg is the (5-function Pg — 5§ s . Introducing 
the Fourier decomposition of delta into Eq. I|56ll one can reproduce the 
projected result (J5SJ). 

The conservation of additive quantum numbers like baryon number, 
strangeness, electric charge or charm is related to the invariance of the 
Hamiltonian under the U(l) Lie group. In many applications it is impor- 
tant to generalize the projection method to symmetries that are related 
with a non-Abelian Lie group G. An example is the special unitary group 
SU(N) that plays an essential role in the theory of strong interactions. Gen- 
eralization of the projection method would require to specify the projection 
operator or generating function. Consequently, the partition function ob- 
tained with the specific eigenvalues of the Casimir operators that fixes the 
multiplet of the irreducible representation of the symmetry group G could 
be determined. 

To find the generating function for the canonical partition function with 
respect to the symmetry group G, let us introduce the quantity Z(g) via 
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Z(g) = Tr[U(g)e-P H }. (57) 

This expression is a function on the group G with U(g) being the unitary 
representation of the group with g C G. The quantity U(g) can be decom- 
posed into irreducible representations U a (g) 



U{g) = Y,U a {9) (58) 

a 

where a is labelling these representations. From Eq. 157|) and 158|) one has 
Z(g)=J2 Tr M9)e- 06 ] 

a 

= EE^«I U a (<9)e-& | ^ Q ,Q (59) 

where v a labels the states within the representation a and £ Q are degeneracy 
parameters. 

Introducing the unit operator 1 =|)(| into the above equation the ex- 
pression factorizes 

Z(g) = (""'^a I ^(ff) I V a ,i a ){v a ,^ a I e _/3ff | I/ a ,fa) 

= EE("«I Ca(ff) I y a )(ta I e"^ | ta), (60) 

where we have used that, due to the exact symmetry, the only non-vanishing 
matrix elements of e~@ H are those diagonal in v a . The matrix elements of 
U a (g) are only non-zero if they are diagonal in £ Q . Finally, the matrix 
elements of the Hamiltonian are independent of the states within represen- 
tation (since due to symmetry they are dynamically equivalent) and those 
of U (g) of degeneracy factors (since U(g) does not distinguish dynamically 
different states that transform under the same representation). 
The last two sums in Eq. (|60J) can be further simplified as 

£> Q | U a (g) | v a ) = Tr a [U a (g)} = X *(g). (61) 



Particle Production in Heavy Ion Collisions 



39 



The quantity \ a is by definition the character of the irreducible U a (g) 
representation and 

I e- pA I U = wi^ a e~ p& = ^r^(r, V), (62) 

where Z Q is the canonical partition function with respect to the G symme- 
try of the Hamiltonian and d(a) is the dimension of the representation a. 
Calculating Z a one considers under the trace only those states that trans- 
form with respect to a given irreducible representation of the symmetry 
group. 

We have thus connected, through Eq. l|6Ufl and l|61H62f> . the canonical 
partition function with the generating functional on the group 



^)=£^^CT,F) (63) 

The canonical partition function is the coefficient in the cluster decom- 
position of the generating function with respect to the characters of the 
representations. 

The character functions satisfy the orthogonality relation 



1 



d(a) 



dn{g)x* a (9)Xi{9) = $ a ,~< (64) 



where dfi(g) is an invariant Haar measure on the group. 

The orthogonality relation for characters allows to find the coefficients, 
the canonical partition function, in this cluster decomposition. From Eq. (|63[) 
and H64J) one gets 

Z a (T, V) = d(a) J d^g) X * a (9)Z(g) (65) 

This result is a generalization of Eq. I|55|) to an arbitrary symmetry group 
that is a compact Lie group. The formula holds for any dynamical system 
described by the Hamiltonian H . 

To find the canonical partition function we have to determine first the 
generating function Z(g) defined on the symmetry group G. If the symmetry 
group is of rank r, then the character of any irreducible representation are 
the functions of r variables {71,.... , 7r}- Denoting as Jk the commuting 
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generators of G with k = 1, ...,r the character function 

Xa(7i>-!7r) = I e^i=i liJi | t/ Q ) (66) 

is obtained. Here, ^q, labels the state within the representation a. With the 
above form of the characters we can write Eq. (|63[1 as 

Z( 7l ,.., 7r )-Tr[ e ^ + ^:::^ J '] (67) 

Through the Wick rotation tj = —i{3fj,i the generating function Z is just 
the GC partition function with respect to the conservation laws given by 
all commuting generators of the symmetry group G. 

The equations (|65|l and (|67|l are the basis that permits to obtain the 
canonical partition function for systems restricted to any symmetry. The 
simplicity of the projection formula l|65|l is that the operators that appear 
in the generating function are additive they are generators of the maximal 
Abelian subgroup of G. Thus, the problem of extracting the canonical par- 
tition function with respect to an arbitrary semi-simple compact Lie group 
G is reduced to the projection onto a maximal Abelian subgroup of G. 

The calculation of the generating function from the Eq. 1)67(1 can be 
done applying standard perturbative diagrammatic methods or a mean field 
approach. However, if interactions can be omitted or effectively described 
by a modification of the particle dispersion relations by impl emen ting; an 
effective mass, then the trace in Eq. (|H7)l can be worked ou exactly, 
leading to 

%) = ex p[£^N ( 68 ) 

where 7 = (71,. ., 7 r ) and Z\ = J (gVdp/2n 2 )p 2 exp {—\Jp 2 + rri^/T) is just 
the thermal particle phase-space in Boltzmann approximation belonging to 
a given irreducible multiplet of a symmetry group G. The sum is taken over 
all particle representations that are constituents of the thermodynamical 
system. 

3.2.1. Canonical models with a non- Abelian symmetry 

To illustrate how the projection method described above works, we discuss 
a statistical model that accounts for the canonical conservation of non- 
Abelian charges related with the SU c (N)xUb(1) symmetry with N = 3 
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and B being the baryon number and c denoting the global gauge colour 
symmetry. 

Let us consider a thermal fireball that is composed of quarks and glu- 
ons at temperature T and volume V . We describe the canonical partition 
function that is projected on the global color singlet and exact value of 
the baryon number. The interactions between quarks and gluons are imple- 
mented effectively, resulting in dynamical particle masses that are temper- 
ature dependent, e.g. through m q i g -\ ~ gT. Since the interactions are only 
trivially modifying the dispersion relations one can still use the free particle 
momentum phase-space. Thus, under this assumption, Eq. (|68ll provides the 
correct description of the generating function. The sum in the exponents 
in H68JI gets contributions from quarks, antiquarks and gluons that trans- 
form under the fundamental (0,1), their conjugate (1,0) and adjoint (1,1) 
representation of the SU C (N) x Ub{^) symmetry group. Thus, 

In Z(T, V, 7, 7B ) = + (69) 

where 7 = (71, ..,7jv_i) are the parameters of the SU C (N) and jb of the 
Ub (1) symmetry group. 

Through an explicit calculation of one-particle partition functions for 
massive quarks and gluons the corresponding generating function is ob- 
tained as 

qo m 2 n VT ™ (-1)"+! 

^ n=0 

[ e^ n ' T XQ {ni) + e~^ n / T x* Q M)] (70) 

where the two terms in the bracket represent the contribution of quarks 
and antiquarks, respectively. The corresponding result for massive gluons 
reads 

hxZ G {T,V,^, lB ) = ^^|^f; J-K 2 (m G /T)[ XG (nj) + x* Q M)} (71) 
da ^tt — f\r 

° 71=0 

where g G , go. and oIq = N,d G = N 2 — 1, are respectively, the quark and 
gluon degeneracy factors and dimensions of the representations. 

Now we can apply this generating function in the projection formula 
(|65[1 to get the canonical partition function. Of particular interest is the 
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color singlet partition function that represents global colour neutrality (phe- 
nomenological confinement) of a quark-gluon plasma droplet. The conju- 
gate character for the SU C (N) singlet representation is particulary simple, 
^.(0,0) _ -y baryon number be treated grand canonically requiring a 

substitution 73 = —i/ig/T in Eq. (|70|l . To find Z one still needs an explicit 
form of the fundamental and adjoint characters and the Haar measure on 
the SU C (N) group. Here we quote their structure for the SU C (3) group. The 
real Lr and the imaginary Lj parts of the character in the fundamental 
(quark) representation are 

L R = cos 71 + cos 72 + cos(7i + 72) 

Lj = sin 71 + sin 72 — sin(7i + 72). 
For the adjoint (gluon) representation 

Xg = 2[cos(7i - 72) + cos(27i + 72) + cos(27 2 + 71) + 1]. (72) 
The invariant Haar measure on the SU C (3) internal symmetry group 



, / \ - . ■ 2/ 2 7l+72x . 2 / 2 72 +7U ( *70\ 

dM(7i,72) = ^2 sm ( 2 ' 2 ' *■ 2 ^ ' 

From Eqs. |EQJ-(£U and (JS^ we write the final result for the SU C (3) 
color singlet partition function that for non-vanishing baryon chemical po- 
tential hb reads 

Z°(fi B ,T,V) = J ^(71,72) exp{cixc + 

c 2 [L R cosh(/3/x B ) + iij sinh(/3)] } (74) 

where the constants c\ and C2 can be extracted from Eqs. (|70II71I) . 

The above partition function shows a complex structure of the inte- 
grand. However, due to its symmetry it is straightforward to show that 
the partition function is real. The thermodynamical properties of this color 
singlet canonical partition function and other thermodynamical observables 
can be studied I I Iby a numerical analysis. 

In finite temperature gauge theory the zero component of the gauge field 
Aq takes on the role of the Lagrange multiplier guaranteeing that all states 
satisfy Gauss law. In Euclidean space one can choose a gauge in such a way 
that Aq{x,t)X v is a constant in space-time. In such a gauge the Wilson 
loop defined as 

1 f 13 
L(x) = -TYPexp[igJ^A (x,T)dT} (75) 
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repre sents the character of the fundamental representation of the SU C (N) 

jinn 

The effective potential of the ,S7 ' (N) spin model for the Wilson loop in 
the above gauge 

coincidejnSEEl 

essentially with the generating function 
given in Eq. 1)74(1. In addition, this generating function could be also related 
to the strong coupling eff ective free energy of the lattice gauge theory with 
a finite chemical potentia l^ ^ 1 Thus, the effective model formulated above 
connects the colored quasi-particle degrees of freedom with the Wilson loop. 



3.2.2. The canonical partition function for Abelian charges 

In this section we show how the projection method described above leads to 
a description of particle yields under the constraints imposed by the Abelian 
£/b(1) symmetry. In this case the formalism is particularly transparent due 
to a simple structure of the symmetry group. 

The U(l) group is of rank one, thus the characters of the representations, 
numbered by the eigenvalues of the conserved charge B, depend only on one 
parameter (p. They are of the exponential type: 

= e <B * (76) 

For the conservation of a few Abelian charges inside the system like stran- 
geness (S), baryon number (B) or electric charge (Q) and charm (C) one 
needs to account for the products of the U(l) symmetries: Us {I) x ^b(I) x 
Uq(1) x UciX)- In this case the characters are numbered by the values 
of all conserved charges and they are expressed as the products of the 
corresponding characters of U(l) groups. For simultaneous conservation of 
baryon number and strangeness the characters read: 

x Us(i)xu B (i) ~ Xu s (i) Xu B (i)- e >■ y") 

The invariant measure on the Us(l) x Ub(X) g r °up is just the product 
of the differentials d^(0s, 0b) = {d(f>s /2tt) ■ (d0s/27r). 

In nucleus-nucleus collisions the absolute values of the baryon number, 
electric charge and strangeness are fixed by the initial conditions. Mod- 
elling the particle production using statistical thermodynamics, in general, 
requires a canonical formulation of all these quantum numbers. We restrict 
our discussion only to the case when at most two conserved charges can be 
simultaneously canonical (e.g. the strangeness and baryon number) and all 
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others are treated using the GC formulation. The corresponding canonical 
partition functions can be obtained from Eqs. I|t)5l(i8l) as: 



Zs = h fj d ^ lS ^(T, V, 4>) (78) 

and 

Z B,s = ^j dcfre-wj d^e- iS ^Z{T,V,4>,^) (79) 

where Z is obtained from the grand canonical (GC) partition function re- 
placing the fugacity parameter Xb, X$ by the factors e 1 ^ and e 1 ^ respec- 
tively, 

Z(T, V, 4>) = Z GC (T, V, X B -» e 4 * A s -> e^) (80) 

The particular form of the generating function Z in the above equation is 
model dependent. In applications of the above statistical partition function 
to the description of particle production in heavy ion and hadron-hadron 
collisions we calculate Z in the hadron resonance gas model. In our anal- 
ysis we neglect interactions between a hadron and resonances as well as 
any medium effects on particle properties. In general, however, already in 
the low-density limit, the modifications o f the resona nce width or particle 
dispersion relation could be of importanc e ! ^ 1 9 1 5 H 1 1 5 1 F or the sake of sim- 
plicity, we use a classical statistics, i.e. we assume a temperature and density 
regime such that all particles can be treated using Boltzmann statistics. 

Within the approximations described above and neglecting the contri- 
butions of multi-strange baryons, the generating function in Eq. Ij78|l . has 
the following form 

Z(T,V ) /i Q ,/i B) ^)=exp(iV s= o + iV s=1 e^+iV s= _ 1 e-^) (81) 

where iV s= o,±i is defined as the sum over all particles and resonances having 
strangeness 0, ±1, 

AUo, ±1 = Z k ( 82 ) 

k 

and Zi is the one-particle partition function defined as 
ri _ Vg k 



m k T K 2 {mk/T) exp(S fc ^ B + QkfJ-Q) (83) 
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with the mass mu, spin-isospin degeneracy factor Qk, particle baryon num- 
ber Bk and electric charge Qk- The volume of the system is V and the 
chemical potentials related to the charge and baryon number are deter- 
mined by hq and fis, respectively. 

With the particular form of the generating function (|81|) the canonical 
partition function Z$ is obtained from Eqs. I|78H80|1 as 

Z S = Z ±J^d4>e- iS *e 3 " t *+ s -"**„ (84) 

where Zq = exp (Ns=o) is the partition function of all particles having zero 
strangeness and where we introduce S±i = N s= ±i with N s= ±i defined as 
in Eq. (JHH- 

To calculate the canonical partition function (|84() one c an ex pand each 
term in the power series and then perform the <f> integratio n^® ! Rewriting 
the above equation as 



1 r 27T v /S 1 S- 1 (Jj^e i 1'+, ^-e- i f) 

Zs = Z Q ^ d0e- s *e V V^" > (85) 

27T Jo 

and using the following relation for the modified Bessel functions Is(x), 

eHt+- t )=Y / t S Is(x), (86) 

— oo 

one gets after the (^-integration the canonical partition function for a gas 
with the net strangeness S: 

Z s (T,V^ B , f i Q ) = Z (T,V^B, f i Q )(-^) s / 2 Is(x) (87) 

where the argument of the Bessel function 

x = 2 y /S 1 S- 1 . (88) 

The calculation of the particle density ilk of species k in the canonical 
formulation is straightforward. It amounts to the replacement 

Z\ \ > X k Z\ (89) 

of the corresponding one-particle partition function in equation (|81|) and 
taking the derivative of the canonical partition function (|84|l with respect 
to the particle fugacity Xk 
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d 

n k = TTT lnZs(Afe) 



(90) 



As an example, we quote the canonical result for the density of kaons 
K + and anti- kaons K~ in an environment with a net overall strangeness 
S, 

C Z K+ S -i Is-i(x) n _ Z K- S i Is+l(x) ( s 
HK+ - V Ia(x) ^ - V y/S^Sli Jb{x) ' (91) 



where x = S\ S~y and Z 1 are as in 1(83(1 and (|84|l . 

For the particular case when S\ = S_i the above equation coincide 
with 1)45(1. Thus, the master equation ((42(1 represents the rate for the time 
evolution of the probabilities for which the equilibrium limit corresponds 
to the canonical ensemble. 

The partition function (|85|) and the corresponding results for parti- 
cle densities 1(91(1 were derived neglecting the contribution of multistrange 
baryons to the generating functional 1(81(1 . Multistrange baryons are, how- 
ever, an important characteristics of the collision fireball created in heavy 
ion collisions. Thus, the canonical formalism described above should be ex- 
tended to account for these particles. Under the constraints of the global 
strangeness neutrality condition S — and including hadrons with stran- 
geness co ntent s = ±1, ±2, ±3 the canonical partition function in Eq. ((84(1 
is replacecPCiniby 



4=o = ^ / d4> exp ( £ Sn^A , (92) 

where S n = Y] k Zl and the sum is over all particles and resonances that 
carry strangeness n with Z\ defined as in Eq. 1(83|) . 

The integral representation of the partition function in Eq. 1(92(1 is not 
convenient for a numerical analysis as the integrant is a strongly oscillat- 
ing function. The partition function, however, after <fi integration, can be 
obtained in a form that is free from oscillating terms. Indeed, rewriting 
Eq. lE3 to 

Z C S=0 = ^e s » £ d<f> J] exp (a n e m * + a^e^)] , (93) 
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and using the relation (|86l) one finds, after integratiorJSSI 



oo oo 

Z$ =0 = e s <> E af4«r 2n ~ 3 ^n(x 2 )/ p (x 3 )/- 2n - 3p (x 1 ), (94) 

n— — oo p— — oo 

where 

a,i = V Si/S~i > ^ = 2y/SiS-i (95) 

and /„ are the modified Bessel functions. 

The expression for the particle density, rii , can be obtained from Eq. (|90|) 
and Eq. I|92|) . For a particle i having strangeness s 

Hi = E E a l a 2 a l 2n '' iP ' S I n{x2)Ip{xz)I^2n^p-s{xi). (96) 

■S 1 — Tl— — oo p— — oo 

In the limit of X2 — > and x 3 — > it is sufficient to take only terms with 
n = and p = in Eq. 194fl and l|9t)f, ^Sl In this case the density of particle 
n^ 7 and antiparticle with strangeness content s and s — —s respectively, 
reads 



n^^i^y^ , nf,!^)'^, (97) 
K 6_i J (x) V 6i i (x) 



with x = ^/SiS-i and as in (S3). 

The above equation is an approximation and can be only used for a 
qualitative discussion. The quantitative description of multistrange particle 
production requires the exact result given in Eq. I)94|l and (35) ■ 



3.2.3. The equivalence of the canonical formalism in the grand 
canonical limit 

Discussing the strangeness kinetics in Section 4.1 we have already indi- 
cated that the canonical description of the conservation laws is valid over 
the whole parameter range. The grand canonical formulation, on the other 
hand, is the asymptotic realization of the exact canonical approach. This 
can be indeed verified when directly comparing particle densities obtained 
in the C and GC ensemble. Consider first a thermal system that contains 
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only strangeness 1 particles and their antiparticles. In such an environment 
the GC result for the strangeness s = ±1 hadrons is obtained from l|15|) as 

n^ ±l = %±iA^. (98) 

with the fugacity X s = exp (/ig/T). 

Comparing the above GC and C result of Eq. H91|) with S — one sees 
that 



n c _ GC 

n s=±l — 1 



-a (X) ■ (99) 
where the effective fugacity parameter 

In the limit of large x — > oo the canonical and the grand canonical 
formulations are equivalent. In the opposite limit, however, the differences 
between these two descriptions are large. This can be seen in the most 
transparent way, when directly comparing the two limiting situations of 
the large and small x in the Eq. 1)91(1 . For x — > oo 

lim £M _> i (ioi) 



and the ratio S-%/ -JSiS-i corresponds exactly to the fugacity Xs in the 
GC formulation (|98|l . Indeed, a strangeness neutrality condition in the GC 
ensemble requires that (S) — 0, thus through Eqs. I|80H81|) one has: 



\ S S, - \- L S^ = 0, (102) 

that is X s = S-i/ y/S\ S-\. 

Thus, neglecting multistrange baryons in the generating functional (|92() 
one gets 



Comparing Eq. I|97|) with the GC result l|98|l for the density n s of multi- 
strange particles one finds that 



n c^ n ocl^l (104) 
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However, one needs to remember that the above relation is only valid if a 
thermal phase space of all multistrange hadrons is negligibly small. This 
assumption is, however, questionable particularly when approaching the 
thcrmodynamical limit. 
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Fig. 11. The canonical strangeness suppression factor (see text for the explanation). 
The SPS and the AGS values are shown for Pb-Pb and Au-Au collisions, respectively. 



From Eq. ((103(1 one concludes that the relevant parameter Fs that de- 
scribes deviations of particle multiplicities from their grand canonical value 
reads 

Fs = J M (105) 

The largest differences appear in the limit of a small x where 

lim ^44 -> x/2 (106) 

This limit reproduces the solution of our kinetic equation 1)27(1 for large 
particle number fluctuations. 

The argument of the Bessel functions in Eq. 1(106(1 describes the size of 
the thermal phase-space that is available for strange particles. For a system 
free of multistrange hadrons the argument x can be also identified as being 
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proportional to the total number of strange particle-antiparticle pairs in 
the GC limit. 

The canonical suppression factor Fs(x) is quantified in Fig. Typ- 
ical values of x expected for the SIS, AGS and SPS energies for central 
collisions of different nucleus are also indicated in this figure. Fig. 
shows the importance of the canonical suppression of particle phase-space 
at SIS energies where it can even exceed an order of magnitude. In central 
heavy ion collisions at the AGS and particularly at higher energies (SPS, 
RHIC, LHC), the canonical suppression is seen in Fig. (|llfl to be negligible. 
Thus, here the GC formalism is adequate. In general, the canonical statis- 
tical interpretation of the particle production in central heavy ion collisions 
is important if the CMS collisions energy per nucleon pair becomes less 
than about (2 — 4) GeV. However, as we show later, the canonical suppres- 
sion effect can also be important at high energy for non- central heavy ion 
collisions or for the description of heavy quark production. 

At the end of this section we also formulate thermodynamics of the 
canonical ensemble with strangeness and baryon number being exactly 
conserved. The corresponding partition function was already presented in 
Eq. Q79J1 . Neglecting the contribution from multistrange baryons and the 
particle-antiparticle charge asymmetry, the generating functional Z(T, V, <fi, iff) 
in Eq. (|8*0f reads 



+ 27V s=0: f,=i cos <f) + 2iV s= i ib= i cos(cj) - ip) 

where N s .b is defined as the sum over all particles and resonances having 
the strangeness s and baryon number 6, 



and Z\ — m\ T K 2 (mf : /T) is the thermal phase-space available for a 
particle that carries strangeness s and baryon number b. 

The above generating functional can be applied in Eq. (|79|l to get the 
canonical partition function of the fireball with net value of strangeness S 
and baryon number B 



In Z S (T, V, n Q = 0, <j>, ip) = iV s= o,&=o + 2iV s= i i6=0 cos ip 



(107) 




(108) 



k 



Zb,s(T,V) 



Zq{T,V) 



L 



d<fie 



exp [zn cos 4>] 



(109) 



4tt 2 
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/ dipe~ lS ^ ex.p[zK cos-0 + zy cos(0 — ip)]d(j)dip 
Jo 

where Z = cxp(iV s=0 ,6=o), z K = N s= i ib=Q , z N = N s=0 ,b=i and z Y = 
^V s =i,6=i< This notation indicates the type of particles that carry corre- 
sponding quantum numbers: charge neutral hadrons, strange mesons, non- 
strange and strange baryons. In the exponent of the ip integral we write 

Zk cos Ip + Zy cos(0 — tp) = z(<j>) cos(V> — o>(cp)) (HO) 

where a(—<j)) = —a((f>) and 

z{4>) = {z 2 K + 2z K z Y cos 4> + Zy) 1 / 2 

e«*(*) = ^ + -^e* (Ill) 

Since the ip integral goes over the whole period, we may shift the integration 
by a and perform the tp integral exactly to yield 

Z BS (T,V) = Z ^ T ' V ^> r cos(B<p + Sa(<p)) exp[2z N cos cp] (112) 

7T 

I s (2z{<t>))d4> 

However, the <p integration cannot be solved analytically. 

Starting from the above partition function one can find the mean mul- 
tiplicity of particle species i. To get it one simply (i) separates for these 
species the particle and anti-particle term in Eq. I|1(J8[I . (ii) multiplies the 
relevant one by A, (iii) differentiates with respect to A following Eq. (|9()(l 
and puts A = 1 afterwards. The result for the particle i with the strangeness 
sSi and baryon number readJUSl 



(N^s. = Zl(m u T,V) B - B - S :*> ' V) (113) 

In view of further applications of these results in heavy ion collisions 
we restrict the discussion only to non-strange systems, that is these with 
overall strangeness 5 = 0. 

The results of Eq. (| 1 1 3|) should coincide with the GC value in the limit 
of large B and V, however, with a fixed baryon density B/V. This can be 
showrP^Sl explicitly using a Chebyshev approximation of the corresponding 
integrals. 
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4. The canonical statistical model and its applications 

The results discussed in the last section indicate that the major differ- 
ence between the C and GC treatment of the conservation laws appears 
through a str ong su ppression of thermal particle phase-space in the canon- 



dependence of particle densities. Deviations from the asymptotic GC limit 
are thus expected to be large for low temperature f and/or small volume. In a 
thermal fireball created in heavy ion collisions these parameters are related 
to the CMS collision energy and the number of participating (wounded) 
nucleons, respectively. It is thus clear that the canonical formulation of 
quantum number conservation should be of importance in low energy cen- 
tral and high energy peripheral heavy ion collisions as well as in hadron- 
hadron collisions. In the following section the applications of the canonical 
statistical model in the above collision scenarios and for different conserved 
quantum numbers will be discussed. A special case is the production of 
hadrons containing charm quarks. This will be dealt with in Section 5.3. 

4.1. Central heavy ion collisions at SIS energies 

The number of strange particles produced in heavy ion collisions depends 
on the energy and centrality of the collision. In low energy A-A collisions 
in the SIS/GSI energy range from 1 to 2 A-GeV, the average number of 
strange particles produced in an event is of the order of 10~ 3 . Thus, follow- 
ing the kinetic analysis presented in Section 4.1, a statistical description 
would require the canonical treatment of strangeness conservation. How- 
ever, the conservation of baryon number and isospin can be treated grand 
canonically. Consequently, one expects a different centrality dependence of 
strange and non-strange particle yields. Fig. (|12|l shows experimental data 
on K + and ir + yields divided by the number of participants j4 par t as a 
function of A pax t measured ^ ^ in Au-Au collisions at beam kinetic energy 
of 1 A-GcV. The data indeed exhibit the behavior expected in the canonical 
statistical model: a strong increase of the K + yield per participant and an 
almost constant ir + yield per participant with centrality. 

In the canonical model the particle densities depend on four parame- 
ters: the chemical potentials, /xq and /is, related with the GC description 



ical approach 




as well as in an explicit volume 



The temperature T should be low relative to the lowest particle mass that carries the 
conserved charge. 
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<A part > 

Fig. 12. The yield of kaons and pions measured in Au-Au collisions at beam kinetic 
energy of 1 A GeV from Ref. 11181 vers us ce ntralitv given by the number of participants 
Apart- The line represents the A^ art fij ^^^ lto experimental data 



of the electric charge and baryon number conservation, the temperature T 
and the volume parameter appearing through the canonical treatment of 
the strangeness conservation. Constraints on these variables arise from the 
isospin asymmetry measured by the baryon number divided by twice the 
charge, B/2Q. For an isospin symmetric system this ratio is simply 1, for 
Ni+Ni it is 1.04 while for Au+Au this ratio is 1.25. When considering par- 
ticle multiplicity ratios we are thus left with three independent parameters. 
The volume parameter V that is responsible for the canonical suppression, 
the freeze-out temperature T and freeze-out baryon-chemical potential [Ib 
of the fireball. 

Figs. 116(1 show the location and sensitivity of the freeze-out pa- 
rameters for different particle ratios in the (T — fis) plane when varying 
these ratios in the range obtained at SIS energy. The deuteron to proton 
d/p and the tt + /p ratios provide a good determination of the range of ther- 
mal parameters. The ir + /p curve in the (T — /is) plane shows temperature 
saturation for a large /ib that establishes the upper limit of the freeze- 
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Fig. 14. Lines of constant dcu- 
tron /pr oton ratio in the T—[j,b plane 
obtained^! in the statistical model. 



out temperature T. On the other hand, the d/p ratio fixes the range of the 
freeze-out value for /ib as it shows a steep dependence on the temperature. 6 
The K + /K~ ratio in Fig. I|15fl exhibits a similar behavior as d/p and is 
also independent of the volume parameter as is evident from Eq. I|45|l when 
requiring that S = 0. 

The variation of thermal parameters with the system size is shown in 
Fig. H16J) using as an example the K + /ir + ratio. For a given system size the 
K + /it + ratio clearly determines a lower limit of the freeze-out temperature 
as it saturates for a large fis- Changing the volume parameter V — A/SirR 3 
implies a substantial modification of the line in the (T — [Ib) plane calcu- 
lated for a fixed value of the K + /ir + yields. Thus, in the canonical model, 
the strange to non-strange particle ratio requires an additional consider- 
ation of the range of correlation of strange particles that is quantified by 
the volume parameter V. This parameter is assumed to be related to the 
number of nucleons participating in A-A collisions. From a detailed analy- 



s The deuteron, is the composite object as it is the proton-neutron bound state. It is most 
likely produced by nucleon coalescence at kinetic freeze-out. Thus, one could question 
if deuteron yield can be used to fix chemical freeze-out parameters (see Section 3.04). 
At SIS energies, however, chemical and thermal freeze— out coincide and the deuteron 
multiplicity seems to follow a statistical order with the same thermal parameters as its 
constituents. 
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sis of experimental data from SIS up to AGS energies it was shown that V 
can be identified^ as the initial overlap volume of the system created in 
A-A collisions. Thus, it is obtained from the atomic number of the colliding 
nuclei and from the impact parameter by simple geometric arguments. In 
heavy ion collisions at SIS energies a good description of Ni-Ni and Au-Au 
data was obtainecP^ h when choosing V ~ VoA part /2 with Vq ~ 7fm 3 i.e. 
of the same order as the volume of the nucleon. 

The comparison of the thermal model with experimental data from 
AGS, SPS up to RHIC energy was discussed in Section 3 and it was shown 
that there are common freeze-out parameters which describe simultane- 
ously all measured particle multiplicity ratios. In order to illustrate that 
this is also the case in low energy heavy ion collisions we show in Fig. I|17|l 
the lines in the (T — Up ) plane corresponding to different particle multiplic- 
ity ratios measured '*"'^ in Ni-Ni collisions at 1.8 AGeV. The experimental 
errors are for simplicity not shown in the figure. All lines, except the one 



h This volume parameter V can be in general y/s dependent. One way to include this de- 
pendence would be to replace the spherical symmetric V by cylinder with its longitudinal 
size being Lorentz contracted. 
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Fig. 18. A com parison of the statistical 
model resultj^21 for the K^/Apart ratio 
with the data from Fig. 1121 . The dashed 
and dashed-dotted lines represent predic- 
tions of the statistical model without and 
with a small A par t dependence of fig and 
T. For more details see text. 



for 77/71"°, have a common crossing point around T ~ 70 McV and \xb ~ 760 
MeV. A value of R ~ 4 fm is needed to describe the measured K + /tt + ratio 
with the freeze-out parameters extracted from ir + /p and d/p ratios. This 
radius is compatible with that expected for a central Ni-Ni collision and 
was founcPSlto be the same in the whole energy range from 0.8 up to 1.8 
A- GeV. 

The corresponding results for a thermal description of Au-Au collisions 
at two different incident kinetic energies 1.0 and 1.5 AGeV can be found 
in Ref. I|49|l. As for Ni-Ni data, the particle ratios, tt + /p, K + /ir + , ir + /ir~ , 
K + /K~ and d/p, with exception of t]/tt , could be described with the same 
value of the freeze-out parameters. The temperature T ~ 53 MeV and 
fj,B ~ 822 MeV were found in Au-Au collisions at 1.0 A-GeV.' Thus, the 
freeze-out temperature is obviously decreasing whereas the baryon chemical 
potential is an increasing function of the collision energy. A small variation 



'With these thermal parameters the total density of particles at chemical freezeout cor- 
responds to ng ~ 0.05/fm 3 . Due to rather large value of the baryochemical potential 
and leading contribution of baryons to total particle density this value also corresponds 
to the total baryon density of the system. 
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of freeze-out parameters with A in central A-A collisions for the same 
collisions energy was also extracted from the datcl^l 

The observed scaling of the volume parameter that determines the 
canonical suppression of the strange particle phase-space with the num- 
ber of participants was also found to be valid in Au-Au collisions. Here, in 
the most central collisions a radius of ~ 6.2 fm is required to reproduce the 
measured K + /A part and K~ /A part yields. The larger radius value obtained 
for Au compared to Ni data is compatible with the larger size of Au and 
corresponds to A part ~ 330 in the most central Au-Au collisions. 

The importance of strangeness suppression due to canonical treatment 
of the conservation laws is particularly transparent in the comparison of the 
thermal model with the Au-Au data at 1 A-GeV . Using the grand canon- 
ical formulation of the strangeness conservation one would get a value of 
K + /n + ~ 0.04 that overestimates the data by more than an order of mag- 
nitude. This shows that the thermal particle phase-space at SIS energies is 
far from the grand canonical limit and that the exact and local treatment 
of strangeness conservation is of crucial importance. 

The multiplicity of K + per participant (K + /A part ) was indicated in 
Fig. 112(1 to increase strongly with centr a lity w hile the corresponding ratio 
for pions tt + /A partl is constant with A par P^. Consequently, the pion yield 
is proportional to the number of participants while the multiplicity of K + 
scales with A part as K + ~ A part with a ~ 1.8. The canonical treatment 
of the strangeness conservation predicts the yield of strange particles to 
increase quadratically with the number of participants, see e.g. Eq. I|28|l. 
In Fig. 1(18(1 the experimental data from Fig. 1)12(1 are compared with the 
thermal model. The parameters were chosen to reproduce the ir/p and d/p 
ratios. The dashed-line in Fig. 1(181) describes the results of a thermal model 
under the assumption that both T and \ib are independent of A part . One 
sees that already under this approximation the agreement of the model 
and the experimental data is very satisfactory as the model describes the 
magnitude and the centrality dependence of these data. Some differences 
between the model and the data seen in Fig. (|T%)) can be accounted for 
when including a small variation of the freeze -out parameters with A part . 
A smooth and almost linear increase of the temperature with centrality by 
a few MeV and a corresponding decrease of fiB from very peripheral to 
central collisions is sufficient to get a very good description of the data. In 
Fig. ((18(1 the dashed-dotted line shows thermal model results that include 
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the variation of thermal parameters with centrality. Such a small depen- 
dence of the freeze-out temperature on impact parameter comes at first 
glance as a surprise since the experimental result on the apparent inverse 
slope parameter T app of particle yields a s a funct ion of pt in Au-Au colli- 
sions shows a strong dependence on -Ap ar J^ ^^ I This difference, however, 
can be accounted for by including the concept of centrality dependent trans- 
verse collective flow of the collision fireball. A detailed analysis of particle 
spectra in Au-Au collisions at 1 A-GeV has shown that keeping the chem- 
ical freeze-out temperature at 53 MeV and including collective transverse 
flow reproduces the transverse momentum distributions of pions, protons 
and kaons as well as their centrality dependencd^. This result indicates 
that chemical and thermal freeze-out coincide at SIS energy. 

The canonical thermal model provides a consistent description of the ex- 
perimental data in the GSI/SIS energy range. The abundances of K + , K~ , 
p, g?, 7r + and 7r~ hadrons (with the notable exception of r\ and possibly <fi) 
J seem to come from a common hot source and with well defined tempera- 
ture, T w 50, 54, 70 MeV, and baryon chemical potential \i B ~ 825, 805, 750 
MeV for central Ni-Ni collisions at 0.8, 1.0, 1.8 A-GeV and correspondingly 
T w 53 MeV and [i B w 822 MeV for central Au-Au collisions at 1.0 A-GeV. 
These temperatures are lower than the ones observed in the particle spec- 
tra but here a common explanation is possible in terms of hydrodynamic 
flow. The flow differentiates between particles of different mass since they 
acquire the same boost in velocity but very different boosts in momentum. 

The common freeze-out condition for almost all particle species is very 
strong evidence for chemical equilibrium in low energy heavy ion collisions. 
A satisfactory agreement of the model and the data could be only obtained 
when including a canonical description of the strangeness conservation. The 
relevant parameter that quantifies the canonical suppression of the particle 
phase-space was found to be the initial volume of the collision fireball that 
scales with the number of projectile participants. This volume parameter 
describes the range of strange particle correlations and is smaller than the 



■■The observed deviations of the model from the measured yield of r) mesons requir»I=2] 
further studies. It is conceivable that the canonical model described in Section 4 does not 
account correctly for hidden st range particle production. The recent result of (j>/K~ = 
0.44±0.16±0.22 ratio obtainecPIZl in Ni-Ni collisions at 1.8 A-GeV could also be larger 
then the statistical model value of 4>IK~ ~ 0.1. However, within large experimental 
uncertainties the model prediction is still not excluded. 
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radius of the fireball at chemical freeze-out. In central Au-Au collisions at 
1 A-GeV the correlation radius was found to be 6.2 fm roughly correspond- 
ing to the size of Au whereas the radius req uired to reproduce measured 
particle yields is almost two times larger^. The appearance of two dif- 
ferent space-like scales in the canonical description of particle production 
can be possibly understood from the kinetics of strangeness production. 
Introducing the locality of strangeness conservation in the kinetic equation 
H15fl implies that the volume parameter in the loss terms can be different 
(smaller) than in the gain terms. Consequently, strange particle yields de- 
pend on two volume parameters: (i) the volume of the fireball which is also 
an overall normalization factor that determines the total strange and non- 
strange particle yields originating from the collisions fireball and (ii) the 
volume that parameterizes the space-like correlations of strange particles 
that is required to satisfy exactly the strangeness conservation. The second 
parameter is also related with the initial number of nucleons participating 
in the collision. 

The apparent chemical equilibration of particle yields measured at SIS 
energies and the kinetic theory developed in Section 4 has recently inspired 
a more complete dynamica l study of the problem in terms of microscopic 
transp ort mode l J^^ W^ffi Recently the relativistic transport model was 
applieeff ^l^ to describe the chemical equilibration of kaon and anti-kaon 
at energies that are below the N-N threshold. The results of these micro- 
scopic studies indicate that K + can possibly appear k in chemical equilib- 
rium during the lifetime of the collision fireball. The K~ , on the other hand, 
approaches chemical equilibrium even at a earlier times However, it may 
eventually fall out of equilibriu m at a later time due to the large annihilation 
cross sections in nuclear matteffini Thus, the results of transport models 
do not exclude chemical equilibration in low energy heavy ion collisions. 
The level of equilibration in these models is strongly related to the magni- 
tude of prod uction and absorption cross section of kaons inside the nuclear 
medium I Although significant progress has been made in the theoret- 
ical description and understanding of in-medium kaon cross sections, the 
results are still far from complete. In particular, recently it was suggested 
that the coupling of kaons with the p-wave £(1385) resonance can substan- 



k However, this equilibration can be obtained if the if + mass is substantially changed in 
the hot and dense medium. 
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tially increase the K production cross section that could influenc e ' the 
approach of kaons towards an equilibrium. 

4.2. Particle production in high energy p—p collisions 

The success of the statistical approach in the description of particle pro- 
duction in heavy ion collisions discussed in Sections 3 and 4 prompts the 
question if a statistical order of the secondaries is also observed in elemen- 
tary collisions such as high energy p-p interactions, to which we restrict 
our discussion. The results of Section 4 make it clear that one should in 
this case apply a model that accounts for the canonical conservation of the 
quantum numbers. This is particularly the case for strangeness since even 
in very high energy p-p collisions the number of produced strange particles 
per event is of the order of unity. Thus, large, event by event, strange par- 
ticle multiplicity fluctuations prevent the applicability of the GC approach. 
Whether or not the GC treatment of the baryon number and electric charge 
is adequate, would require a detailed study of the relative error between the 
canonical and the grand canonical results. 

The application of the statistical model to the eleme ntary hadron- 
hadron reactions was first proposed by Rolf HagedorrPEI i n order to de- 
scribe the exponential shape of the mt-spectra of produced particles in p-p 
collisions. Hagedorn also pointed out phenomenologically the importance 
of the canonical treatment of the conservation laws for rarely produced 
particles. The first application of the canonical model to s trangeness pro- 
duction in p-p collisions was done by Edward ShuryalfiSSI i n the context 
of ISR data. Recently a complete analysis of hadron yields in p-p as well 
as in p-p, e + e _ , 7r-p and in K-p collisions at several center-of-mass ener- 
gies has been done in Refs. I|451 I56lll33[l . This detailed analysis has shown 
that particle abundances in elementary collisions can be also described by 
a statistical ensemble with maximized entropy. In fact, measured yields are 
consistent with the model assuming the existence of equilibrated fireballs 
at a temperature T «160-180 MeV. 

The most general partition function Zb,q,s{Vi T) that is applied to test 
the chemical composition of the secondaries in elementary collisions should 
account for the canonical conservation of baryon number B, strangeness S, 
and electric charge Q. It can be constructed applying the projection method 
(see Section 4.2) for the U B (l) x Us (I) x Uq(1) symmetry. Following Eq. (JEJ 
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Fig. 19. A comparisor P^ I of the p-p multiplicity data at y/s = 27 .4 with the statis- 
tical model that accounts for the exact conservation of baryon number, electric charge, 
strangeness and includes the strangeness undersaturation factor 7 S ~ 0.51. 

one has 



(V,T) 



2n 



2tt 

d(j>s 
2tt 



2tt 



-iQ4>Q 



e -iS^ Se \nZ(Ty,4, B ,<t>Q,<t>s) 



(114) 



where B,Q, and S are the initial values of the quantum numbers in hadron- 
hadron collisions which are in p-p scattering B=Q=2 and S=0. 

The generating function Z in Eq. I|114|) can be obtained from the grand 
canonical partition function of a hadron resonance gas © by a Wick rota- 
tion of all appropriate chemical potentials 



Z(T,V,4>b,4>q 
Z GC (T,V,fi B 







Ms -> -#s). (H5) 



Different particle multiplicities and their ratios are obtained from the 
partition function (|114(l following the procedure that was described in Sec- 
tion 4.2.3. 
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Fig. <|19|l shows an example of the comparison of the canonical model 
prediction 1114|) with the experimental data for different particle yields 
obtainecP>22l hi p-p collisions at i/s — 27.4-27.6 GeV. The agreement of 
the model with most of the experimental data is reasonable. However, the 
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Fig. 20. A comparison of the p-p mul- 
tiplicity obtained at = 27.4 with the 
canonical models. The open symbols rep- 
resent the result of the model that accounts 
for S, B and Q being exactly conserve 
The filled symbols are the statistical model 
results with only S being exact whereas B 
and Q are treated grand canonically. 



Fig. 21. As in Figure |2UJ but the filled 
circle are obtained in the canonical model 
that accounts for the strangeness undersat- 
uration to be controlled by the correlation 
volume, V p = 4nR p /3 with Rp ~ 1.1 fm 
instead of 7 S . 



yields of the resonances like A°,p° and 4> can differ by a few standard 
deviations from the data leading to rather large reduce \ 2 values of ~ 5. 

The agreement of the model with the strange pa rticle yiel ds could only 
be, however, achieved in Ref. when introducing22H22l an additional 

parameter j s into the canonical partition function l|114Jl . This parameter 
suppresses a thermal phase-space of particles composed of n s strange or 
antistrange quarks by a factor (j s ) ns . In p-p and p-p collisions and in a 
very broad energy range from yfs ~ 20 up to yfs ~ 900 GeV the same value 
of 7 S ~ 0.51 was needecP^ to reproduce the measured strange particle 
yields. 

In high energy p-p collisions the canonical effects due to the baryon 
number and isospin conservation are expected to be small. Fig. l(2T)jl shows 
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the comparison of p-p data with the canonical model that accounts for 
only strangeness being conserved exactly. The baryon number and electric 
charge conservation are treated in the GC ensemble, thus are controlled by 
chemical potentials. A satisfactory description of experimental data with 
Hb ~ 195 MeV, fi Q ~ 30 MeV, T ~ 165 MeV and 7s - 0.53 seen in 
Fig. H20|) shows that canonical e ffects related with charge and baryon num- 
ber conservation are indeed smalff^]. This is, however, not the case in p-p 
and e + e~ collisions since the initial values of B = Q = are there obviously 
too small to use a GC approximation. 

The experimental data shown in Fig. (|19|l can be also described in terms 
of the canonical model that was successfully applied in low energy heavy ion 
collisions (see Section 5.1). There, instead of the strangeness undersatura- 
tion factor 7,5, space- like correlations of strange particles were introduced. 
Consequently, there were two volume parameters that determined parti- 
cle yields: (i) the volume of the fireball and, (ii) the correlation volume. 
Fig. H21|l shows that choosing the correlation volume V p — and 
Rp ~ 1.1 — 1.2 fm the p-p data are well reproduced with the exception of 
<fi. The </>-meson is not canonically suppressed because it only carries zero 
strangeness. Taking the correlation volume instead of "f s gives a deviation of 
the measured (f> abundance from a thermal fit of 12a, thus increasing even 
more the already large deviation of 4 cr se en in Fig. (|19|1 . The canonical 
model that introduces a j s ~ 0.5 factor^ 

is obviously better reproducing 
the yields of hidden strange particles obtained in elementary collisio ns. 

The value of 7 S extracted from p-p data is smaller than required^ in 
central Pb-Pb at 160 A-GeV where for fully integrated yields 7 S = 0.75 ± 
0.05 was fitted 1 . This result shows that the undersaturation of strangeness in 
p-p and A-A collisions differs by almost 50%. This difference alone already 
indicated^ that strangeness in A-A is enhanced relative to p-p collisions. 111 

The temperature parameter extracted from particle yields in high en- 
ergy elementary collisions is at first glance surprisingly compatible with the 



'We have to point out, however, that midrapidity data in A— A collisions at SPS are 
consistent*''! with the value of 7 S = 1. See the more detailed discussion in Section 2.2. 
m Looked at from a different angle we conclude that in central nucleus— nucleus colli- 
sions at AGS energy and higher strangeness production reaches values consistent with 
complete chemical equilibrium. At lower energies, and in particular in elementary par- 
ticle collisions, strangeness is strongly undersaturated. This will be discussed further in 
Section 6. 
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chemical freeze-out temperature extracted from heavy ion datcPnEHl At 
SPS and RHIC energies, T can be considered as a measure of thermal exci- 
tations of a non-perturbative QCD vacuum due to the particle scattering. 
Thus, T should be mostly correlated with the collision energy and not with 
the system size. The charge chemical potential, however, due to the isospin 
asymmetry, differers substantially in p -p a nd A-A collisions. In Pb-Pb col- 
lisions at the SPS the fiQ ~ —7 Me\*22l whereas in p-p at y/s ~ 27 the 
(1q ~ 35 MeV is required. 



4.2.1. Statistical hadronization and string dynamics in p-p 
collisions 

The apparent agreement of the canonical statistical model with experimen- 
tal data on particle production in elementary collisions leads to the inter- 
pretation that hadronization in particle collisions is a statistical process. 
This result is difficult to reconcile with the popular picture that hadron 
production in hadron hadron collisions is due to the decay of color flux 
tubeJi^l a model that has explained many dynamical features of these 
collisions. In the following we address the question, how one can possibly 
distinguish the string hadronization via the break up of a color flux tube 
from the statistical hadronization. We argue following Ref. I|136[) that the 
il/Q = f2 + /f2~ ratio in elementary proton-proton collisions is a sensitive 
probe to differentiate possibly these two scenarios. 

Color flux tubes, called strings, connect two SU(3) color charges [ 3 ] and 
[ 3 ] with a linear confining potential. If the excitation energy of the string 
is high enough it is allowed to decay via the Schwinger mechanisnfflZI i.e. 
the rate of newly produced quarks is given by: 

— — - ~ exp [— ixm\l k\ , (116) 
dp± 

where k is the string tension and mi = \Jp\+ m 2 is the transverse mass 
of the produced quark with mass m. 

However, specific string models may differ in their philosophy and the 
types of strings that are created: 

• In UrQMEDUD 

the projectile and target protons become excited 
objects due to the momentum transfer in the interaction. The re- 
sulting strings, with at most two strings being formed, are of the 
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diquark- quark type. 

• In NeXia I the p-p interaction is described in terms of pomeron 
exchanges or ladder diagrams. Both hard and soft interactions take 
place in parallel. Energy is equally shared between all cut pomerons 
and the remnants. The endpoints of the cut pomerons (i.e. the 
endpoints of the strings) may be valence quarks, sea quarks or 
antiquarks. 

• In PYTHI^PZI 

a scheme similar to that in UrQMD is employed. 
However, hard interactions may create additional strings from scat- 
tered gluons and sea quarks. Most strings are also of a diquark- 
quark form. 

Fig. |22I -left) depicts the antibaryon to baryon ratio at midrapidity in 
proton-proton interactions at 160 GeV. The results of th e cal culations by 
NeXuS, UrQMD and PYTHIA are included in this figurJEi In all these 
models, the B/B ratio increases strongly with the strangeness content of the 
baryon. For strangeness \s\ = 3 the ratio significantly exceeds the unity. In 
UrQMD and PYTHIA the hadronization of the diquark-quark strings leads 
directly to the overpopulation of tt. In NeXuS, however, the imbalance of 
quarks and anti quarks in the initial state leads to the formation of q va i~ s soa 
strings (the s va i — 9 S ca string is not possible). These strings result then in 
the overpopulation of f2's. 

In order to understand the large f2/fi values predicted by string models 
we include in Fig. I|23[l . the color flux tube break-up mechanism. 
Fig. (|23|l shows the fragmentation of the color field into quark-antiquark 
pairs, which then coalesce into hadrons. While in large strings fi's and O's 
are produced in equal abundance (a) , low-mass strings in UrQMD suppress 
production at the string ends (b), while in NeXuS Q's are enhanced (c). 
Thus, the microscopic method of hadronization leads to a strong imbalance 
in the f2/f2 ratio in low- mass strings. 

In Fig. (|22[1 the string mo del re sults are compared with the predictions 
of two s tatis tical models (SM^P^ anc [ the preliminary experimental data 
obtainecP^by the NA49 Collaboration. The main difference between these 
models is the implementation of baryon number and electric charge conser- 
vation and the way an a ddit ional strangeness suppression is introduced. In 
model (i) the calculation^! is a full canonical one with fixed baryon num- 
ber, strangeness and electric charge identical to those of the initial state. An 
extra strangeness suppression is introduced to reproduce the experimental 
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Fig. 22. The left hand figure: the anti baryon to baryon ratio at \y — y C m\ < 1 in 
p-p interactions at 160 GeV as given by PYTHIA, NeXuS and UrQMD. The right hand 
figure: the anti baryon to baryon ratio for the same reaction as given by the statistical 
models. Stars depict preliminary NA49 data for the B/B ratio at midrapidity. 

multiplicities. This is done by considering the number of newly produced 
(ss) pairs as an additional charge to be found in the final hadrons. The ss 
pairs fluctuate according to a Poisson distribution and their mean number 
is considered as a free parameter to be fitted^. The parameters used for 
the prediction of the f2 + /f2~ ratio (T, the global volume V sum of sin- 
gle cluster volumes and (ss)) have been obtained by a fit to preliminary 
NA49 p-p datcPi yielding T = 183.7 ± 6.7 MeV, VT 3 = 6.49 ± f. 33 and 
(ss) = 0.405 ± 0.026 with a x 2 /dof = 11.7/9. It must be pointed out that 
the f2 + /fi - ratio is actually independent of the (ss) parameter and only 
depends on T and V. 

In model (ii) the conservation of the baryon number and electric charge 
is approximated by using the GC ensemble. The strangeness conservation 
is, however, implemented on the canonical level following the procedure 
that accounts for strong correlations of produced strange particles. In p-p 
collisions the strangeness is assumed not to be distributed in the whole 
volume of the fireball but to be locally strongly correlated. A correlation 
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Fig. 23. Fragmentation of a color field into quarks and hadrons. While in large strings 
Q-s and O's are produced in an equal abundance (a), small di'quark strings suppress Q's 
at the string ends (b), sea-s quarks enhance f2's (c). 

volume parameter Vq — AttRq/3 is introduced, where Rq ~ 1 fm is a typical 
scale of QCD interactions. The temperature T ~ 158 MeV and /is ~ 238 
MeV were taken as obtained^ from the SM analysis of a full phase-space 
Pb-Pb data of NA49 Collaboration. The volume of the fireball V ~ 17 fm 3 
and the charge chemical potential in p-p was then found to reproduce the 
average charge and baryon number in the initial state. 

The predictions of the statistical models are shown in Fig. 1)221 -right). 
In these approaches the B/B ratio is seen to exhibit a significantly weaker 
increase with the strangeness content of the baryon than that expected 
in the string fragmentation models. For comp ariso n, both figures include 
preliminary data on the B/B ratios obtainecP^I a t midrapidity by the 
NA49 Collaboration. Note that the predictions of the statistical models in 
Fig. Q22JI refer to full phase-space particle yields whilst measurements of 
B/B ratios in p-p collisions have been performed at midrapidity, where 
they are expected to be the largest. Therefore, sizeable deviations of the 
model results from the data seen in Fig. (12211 are to be expected. Howe ver, 
admitting the applicability of SM for midrapidity one reproduced^!! the 
experimental data quite well. 
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In macroscopic string models the fi/Q ratio depends in a strongly non- 
linear fashion on the mass (energy) of the fragmenting string. All these 
models predict a strong enhancement of O production at low energies, while 
for large string masses the ratio approaches the value of O/O = 1 (which 
should be reached in the limit of an infinitely long color flux tube). 

Statistical models, on the other hand, are not able to yield a ratio of 
Cl/Q > 1. This can be easily understood in the GC formalism, where the 
B/B ratio is very sensitive to the baryon chemical potential hb- For finite 
baryon densities and including 100% feeding from resonances, the B/B 
ratio will always be < 1 and only in the limit of /is = may O/O = 1 be 
approached. These features survive in the canonical framewo rk, where the 
GC fugacities are replaced by the ratios of partition function^^EEEIZI 

From the above discussion and from Fig. 12_!l) it is thus clear that within 
the fragmenting color flux tube models the O /O ratio is significantly above 
the unity. This is in strong contrast to statistical model results, that always 
imply that B/B ratios are below or equal to a unity in proton-proton 
reactions. Since this observable is accessible to NA49 measurements at the 
SPS it can provide an excellent test to distinguish the statistical model 
hadronization scenario from that of a microscopic color-flux tube dynamics. 

We have to point out, however, that the classical string models con- 
sidered above do not account for the so called string junction mechanism 
that allows for diffusion of baryon number towards midrapidity. This mech- 
anism, recently includecPH in Dual Parton Model, was shown to be very 
important for (multi)strange baryon and antibaryon production. Thus, it 
would be of importance to study the energy dependence of B/B ratio in p-p 
collisions in terms of the model that includes this baryon number transport. 

4.3. Heavy quark production 

Charm quarks are heavy (m c 3> T c ) thus, thermal production of charm 
quarks and char med h adrons is strongly suppressed in ultra-relativistic 
heavy ion collision J^^ I The situation has been recently discussecP^E^with 
the conclusion that, compared to direct hard production, thermal produc- 
tion of charm quarks can be neglected at SPS energies and is small even 
at LHC energy. However, these investigations led to a new scenario for the 
production of hadrons which contain charm quarks in which production of 
heavy quarks through hard collisions is combined with a statistical proce- 
dure to produce open and hidden charm hadrons at hadronization. This idea 
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of statistical had ronizatio n of charm quark d * '"I has sparked an intense ac- 
tivity in this Initial interest focusscd on the available SPS data 
on J/ip prod uctio n and their interpretation in terms of a conventional sta- 
tistical mode^^l As we show below, these data can be well described, but 
only assuming a charm cross section which is enhanced compared to pre- 
dictions within the framework of perturbative QCD. However, the largest 
differe nces between results from the statistical coalescence scenario (or a 
similai^^ model) and more conventional models a re e xpected at collider 
energies. For example, in the Satz-Matsui approacfP^, one would expect 
very stron g sup pression compared to direct production of J/ip mesons (up 
to a factor ' I of 20) for central Au-Au collision at RHIC energy. In the 
present approach this suppression is overcome by statistical recombination 
of J/ip mesons from the same or different cc pairs, so that much larger 
yields are expected. We therefore focus in this section on predictions 11 for 
open and hidden charm mesons at RHIC and LHC energy, with emphasis 
on the centrality dependence of rapidity densities. 



4.3.1. Statistical Recombination Model 

In this model it is assumed that all charm quarks are produced in pri- 
mary hard collisions and (thermally) equilibrate in the quark-gluon plasma 
(QGP); in particular, no J/ip is preformed in QGP (complete screening) 
and there is no thermal production of charm quarks. For a description of 
the hadronization of the c and c quarks, i.e. for the determination of the 
relative yields of charmonia, and charmed mesons and baryons, we employ 
the statistical mod el, wi th parameters as determined by the analysis of all 
other hadron yield J '^ 1 The picture we have in mind is that all hadrons 
form within a narrow time range at or close to the phase boundary. All 
charmed hadrons (open and hidden) are f ormed at freeze-out (at SPS and 
beyond, freeze-out is at the phase boundary! l^ f) according to the statistical 
laws. 

Another interesting point concerns the ip / (J/ip) ratio. As is well known, 
this ratio i s, in hadron-proton and p-nucleus collisions, close to 12 %, 
independent! I of collision system, energy, transverse momentum etc.. In 
the thermal model, the ratio is 3.7 %, including feeding of the J/ip from 



"This section is based on work by the authors and A. Andronic and reported in Ref. 11481 
Il49l . 
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heavier charmonium states. A temperature of about 280 MeV would be nec- 
essary to explain the ratio found in pp and p-nucleus collisions in a thermal 
approach. Clearly, 3/ip and ip production in pp and p-nucleus collisions are 
manifestly non-thermal. This was previously realized in Refs. I|1541 115l|) . 
Similar considerations apply for the x states. In fact, feeding from xi to 
J/ip is less than 3 % if the production ratios are thermal. 

The experimental situation concerning the evolution with participant 
number of the ip /(3/ip) ratio in nucleus- nucleus collisions, multiplied with 
the respective branching into muon pairs, is presented in Fig. (|24|l . The 
data are from the NA38/50 collaboratio n 1 1 ^ 5 ! 1 ^ ! 15 ? ! 15 ^ With increasing 
Np ar t the ip /(3/ip) ratio drops first rapidly (away from the value in pp 
collisions) but seems to saturate for high N part values at a level very close 
to the thermal model prediction, both for S+U and Pb-Pb collisions 



> 0.02 

-. 

tL 

CO 

£0.015 
> 

=. 

CO 

0.01 



0.005 



• SU data NA50 
■ PbPb95 
□ PbPb96 



pW 
pd 



average pp and pA 



50 100 150 200 250 300 350 



400 450 

Npar, 



Fig. 24. Comparison of the dependence of the measured ip /(J/V 1 ) ratio on the number 
of participating nucleons with the prediction of the thermal model. The data of NA38 
and NA50 Collaborations are from Refs. <15SlllSrjHT571IT^l . See text and Refs. 1751 17M 
for more details. 



Taking this into account we note that predictions of the model should 
only be trusted from about N par t > 150 on, where also the ip' /(J/tp) ratio 
is close to the thermal value for Pb-Pb data. In our approach, ratios for 
all higher charmonia states including the Xc should approach the thermal 
value from N part > 150 on, implying that for those N paT . t values feeding 
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to J/ip should be small. In this picture, there should thus not be different 
"thresholds" for the disappearance of different charmonia. 

The total number of open charm hadrons expected in a purely thermal 
approach , iV™, is then readjusted to the number of directly produced cc 
pairs, jV*' as (neglecting charmonia): N*t = \gJ*%h{gcN%)/ h{g c N%), 
from which the charm enhancement factor g c is extracted. Here, /„ are 
modified Bessel functions. Note that we use here the canonical approach 
as, depending on beam energy, the number of charm quark pairs maybe 
smaller than 1. The grand-canonical limit will likely only be reached at 
LHC energy. For a detailed study of the transition from the canonical to 
the grand-canonical regime see Section 3.2 and Ref. I)76|l . The yield of a 
given species X is then determined by Nx — g c N x h Ii(g c N*^)/ I (g c N^) 
for open charm mesons and hyperons and Nx — g^^x f° r charmonia (see 
Refs. lf7Kll?ri)l for more details). 

The inputs for the above procedure are: i) the total charged parti- 
cles yields ( or rapidity densities), which are taken from experiments at 
SPS^HHUSni an d RffldBHland extrapolated for LHC; and ii) N^ r , which 
is taken from n ext- to - leading order (NLO) perturbative QCD (pQCD) cal- 
culations for (the yield from MRST HO parton distributions was 
used here) and scaled to AA via the nuclear overlap function. A constant 
temperature of 170 MeV and the baryonic chemical potential fif, according 
to the parameterization /^(MeV)=1270/ (1 + y / sjv r /v/4.3 have been used for 
the calculations^! 



4.3.2. Results 
We first compare 

predictions^!!^ 

of the model to 47r-integrated J/tp 
datJSH at the SPS from NA50 Collaboration replotted as outlined in 
Ref. H76J) . In Fig. (|25(l wepresent the model results for two values of N^ r : 
from NLO calculation:J223 (dashed-line) and scaled up by a factor of 2.8 
(continuous line) . The dashed-dotted line in Fig. (|2*5f is obtained with the 
NLO cross section for charm production scaled-up b y a, f actor 1.6, which 
is the ratio of the open charm cross section estimated ™ by NA50 for p- 
p collisions at 450 GeV/c and the NLO values from Ref. Ijl49|l . For this 
case the Nr, n/r t scaling is not the overlap function, but is taken according to 
measured ™ ' dimuon enhancement as a function of N part . 

The results of Fig. f2*K|l indicate that the observed centrality depen- 
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dence of J/ip for 100 < N part <350 is well reproduced by the statistical 
model using the NLO cross sections for charm production scaled by the 
nuclear overlap function. However, to explain the overall magnitude of the 
data a iV* r increase by a factor of 2.8 compared to NLO calculations is 
needed. The drop of the J/tp yield per participant obser ved in the data for 
Npart >350 ( see Fig. Q25fl') is currently understoocPSOl in terms of energy 
density fluctuations for a given overlap geometry. 

We mention in this context that the observed e nhan cement of the dimuon 
yield at intermediate masses has been interprete d^^ ! by NA50 as possible 
indication for an anomalous increase in th e char m cross section. We note, 
however, that other plausible explanation J ^ 1 ^ @ l exist of the observed en- 
hancement in terms of thermal radiation. 



NA50 data 




50 100 150 200 250 300 350 400 
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Fig. 25. The centrality 
dependence of J/ip production 
at SPS. Model predictions are 
compared to 47 r-integrated 
NA50 dati^HR Two curves 
for the model correspond to 
values of N^l r ~ Q 1 37 from 
NLO calculation^!^ (dashed 
line) and scaled up by a 
factor of 2.8 (continuous 
line). The dashed-dotted 
curve is obtained when con- 
sidering the possible NA50 
Npart— dependent charm 
enhancemen^^S] over their 
extracted !*^ ! p-p cross section 
(see text). 



We turn now to discuss model predictions for collider energies. For com- 
parison we include in this study also results at SPS energy. The input 
parameters for these calculations for central collisions (N part =350) are pre- 
sented in Table |3 Notice that from now on we focus on rapidity densities, 
which are the relevant observables at the colliders. The results are compiled 
in Table |31 for a selection of hadrons with open and hidden charm. All pre- 
dicted yields increase strongly with energy, reflecting the increasing charm 
cross section and the concomitant importance of statistical recombination. 
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Also ratios of open charm hadrons evolve with increasing energy, reflecting 
the corresponding decrease in charm chemical potential. 



Tabic 2. Input parameters for model calculations at SPS, RHIC and LHC. 



V s iVJV (GeV) 


17.3 


200 


5500 


T (MeV) 
Mb (MeV) 


170 

253 


170 
27 


170 
1 


dN ch /dy 
V Ay =i (fm 3 ) 


430 
861 


730 
1663 


2000 
4564 


dN^/dy 
9c 


0.064 
1.86 


1.92 
8.33 


16.8 
23.2 


Source: A. Andronic et al. from Ref. 11481. 






Table 3. Results 


of model calculations at SPS, 


RHIC and LHC for N p , 


irt =350. 


y/SNN (GeV) 


17.3 


200 


5500 


dN B+ /dy 
dN u -/dy 


0.010 
0.016 


0.404 
0.420 


3.56 
3.53 


dN u o/dy 
dN A Jdy 


0.022 
0.014 


0.89 
0.153 


7.8 
1.16 


dNj/^/dy 
dN^/dy 


2.55-10~ 4 
0.95-10" 5 


0.011 

3.97-10~ 4 


0.226 
8.46-10 



Source: A. Andronic et al. from Ref. 11481 . 



Predictions for the centrality dependence of J/tp production are pre- 
sented in Fig. (|26(l . In addition to the dramatic change in magnitude (note 
the scale-up by factors of 10 and 100 for RHIC and SPS energy, respectively) 
the results exhibit a striking change in centrality dependence, reflecting 
the transition from a canonical to a grand-canonical regime (see Ref. I|7fi[) 
for m ore d etails). The preliminary PHENIX results on J/tp production at 
RHICJESj agree, within the still large error bars, with our predictions. A 
stringent test of the present model can only be made when high statistics 
J/tp data are available. Another important issue in this respect is the accu- 
racy of the c harm cross section, for which so far only indirect measurements 
are availabl e ! ^*-' ! In any case, very large suppression factors as predicted, 
e.g., by Ref. I|147fl seem not supported by the data. In Fig. I|2t)f) we present 
the predicted centrality dependence of charged D + -meson production for 
the three energies. The expected approximate scaling of the ratio D + / N part 
like Npart (dashed lines in lower Fig. (J2EJ) is only roughly fulfilled due to 



74 



Peter Braun-Munzinger, Krzysztof Redlich, Johanna Stachel 




Fig. 26. Centrality depen- 
dence of J /if) (upper figure) 
and of D+ (lower figure) rapid- 
ity density at SPS, RHIC and 
LHC energies. 



departures of the nuclear overlap function from the simple N p ' art depen- 
dence. 

4.3.3. Charmonium Production from Secondary Collisions at LHC 
Energy 

Another possibility to produce charmonium states is due to reactions among 
D mesons in the hadronic and mixed phase of the collision. This has been 
investigated in Refs. I|142l I167f) . As demonstrated there, this mechanism 
does not lead to appreciable charmonium production at SPS and RHIC 
energies. However, the large number of cc pairs and consequently D,D 
mesons produced in Pb-Pb collisions at LHC energy can lead to an ad- 
ditional production of charmonium bound states due to reactions such as: 
DD* + D*D + D*D* -> J/ip + vr and D*D* + DD — > J/tp + p. These pro- 
cesses were studied within a kinetic model taking into account the space- 
time e volu tion of a longitudinally and transversely expanding medium. The 
result demonstrate that secondary charmonium production appears al- 
most entirely during the mixed phase and is very sensitive to the charmo- 
nium dissociation cross section with co-moving hadrons. Within the most 
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likely scenario for the dissociation cross section of the J/ip mesons their 
regeneration in the hadronic medium will be negligible, even at LHC en- 
ergy. Secondary production of ip' mesons however, due to their large cross 
section above the threshold, can substantially exceed the primary yield. 

4.3.4. Conclusions on Heavy Quark Production 

We have demonstrated that the statistical coalescence approach yields a 
good description of the measured centrality dependence of J/ip production 
at SPS energy, albeit with a charm cross section increased by a factor of 
2.8 compared to current NLO calculation. Rapidity densities for open and 
hidden charm mesons are predicted to increase strongly with energy, with 
striking changes in centrality dependence. First RHIC data on J/tp pro- 
duction support the current predictions, although the errors are too large 
to make firm conclusions. For LHC energies we predict cross sections for 
charm production in central Pb-Pb collisions significantly exceeding the val- 
ues predicted by scaling results for N-N collisions with the nuclear thickness 
function. The statistical coalescence implies travel of charm quarks over sig- 
nificant distances in QGP. If the model predictions will describe consistently 
precision data then 3 ftp enhancement (rather than suppression) would be 
a clear signal for the presence of a deconfined phase. 

Regeneration of charmonia in the mixed and hadronic phase has also 
been studied. For 3/tp mesons this will likely only be a small effect, even at 
LHC energy. However, secondary production of ip' mesons may be signifi- 
cant at LHC energy. 

5. Unified conditions of particle freeze out in heavy ion 
collisions 

A detailed analysis of experimental data in heavy ion collisions from SIS 
through AGS, SPS up to RHIC energy discussed in Section 3 and 5 makes 
it clear that the canonical or grand canonical statistical model reproduces 
most of the measured hadron yields. 

Figure l|27|) is a compilation of chemical freeze-out parameters that are 
required to reproduce the measured particle yields in central A-A collisions 
at SIS, AGS, SPS and RHIC energy. The GSI/SIS results have the lowest 
freeze-out temperature and the highest baryon chemical potential. As the 
beam energy increases a clear shift towards higher T and lower [in occurs. 
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Fig. 27. A compilation of chemical freeze— out parameters appropriate for A-A collisions 
at different energies: SIS results are from Ref. 03B> AGS from Ref. H3i SPS at 40 A-GeV 
from Refs. I44lll68ll4^1ll*75l . SPS at 160 A-GeV from Refs. J2H1EI3II3, and RHIC from 
Refs. 1351 1541 1551 . The full line represents the phenonienological condition of a chemical 
freeze-out at the fixed energ y/particle ~ 1.0 GeV0] 



There is a common feature to all these points, namely that the average 
energy (E) per average number of hadrons (N) is ap prox imately 1 GeV. 
A chemical freeze-out in A-A collisions is thus reached"^ when the energy 
per particle (E)/(N) drops below 1 GeV at all collision energies. 

In cold nuclear matter the (E)/(N) is approximately determined by the 
nucleon mass. For thermally excited nuclear matter, in the non-relativistic 
approximation 



(AO 



(117) 



with (m) being the thermal average mass in the collisions fireball. This 
result makes it clear why at SIS the energy /particle at chemical freeze-out 
is of the order of 1 GeV since T ~ 53 MeV. At SPS and RHIC energy 
the leading particles in the final state (at thermal freeze-out) are pions. 
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However, at chemical freeze-out most of the pions are still hidden in the 
mesonic and baryonic resonances. Thus, here the average thermal mass is 
larger than the pion mass and corresponds approximately to the p-meson 
mass. Consequently, since (m) >> T, Eq. (|117|l can be still used to justify 
approximately that (E)/(N) ~ 1 GeV at the SPS. Actually, this argument 
holds, to a large extent, in the whole energy range from SIS up to RHIC 
energy. 

The physical origin of the phenomenological freeze-out condition of fixed 
energy /particle would require further dynamical justification and interpre- 
tation. Recently, this question has been investigated in central Pb-Pb col- 
lisions at the SPS in ter ms of the Ultra-relativistic Quantum Molecular 
Dynamics model (Ur A detailed study has shown that there 

is a clear correlation between the chemical break-up in terms of inelastic 
scattering rates and the rapid decrease in energy per particle. If (E) / (N) 
approaches the value of 1 GeV the inelastic scattering rates drop substan- 
tially and further evolution is due to elastic and pseudo-elastic collisions 
that preserved the chemical composition of the collision fireball. Following 
the above UrQMD results and the previous suggestions^^ one could con- 
sider the phenomenological chemical freeze-out of (E) / (N) ~ 1 GeV as the 
condition of inelasticity in heavy ion collisions. 

Unified freeze-out conditions were also considere d!^ in the context of 
hydrodynamical models for particle production and evolution in heavy ion 
collisions. There, it was suggested, that the condition for chemical freeze- 
out, (E)/(N) ~ 1 GeV, selects the softest point of the equation of state, 
namely the point where the ratio of the thermodynamical pressure P to the 
energy density e has a minimum. The considerations were essentially based 
on the propose d^ ^ mixed phase model that seems to be consistent with the 
available QCD lattice data. The quantity P/e is closely related to the square 
of the velocity of sound and characterizes the expansion speed EH f the 
reaction zone. Thus, the system lives for the longest time around the softest 
point that allows to reach the chemical equilibrium of its constituents. The 
above interpretation, however, crucially depends on the type of the equation 
of state used in the model. 

Chemical freeze-out in heavy ion collisions can also be determinec P^^ 
by the condition of fixed density of the total number of baryons plus an- 
tibaryons. As it is seen in Fig. (|28l) . within statistical uncertainties on the 
freeze-out parameters the above condition provides a good description of 
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Fig. 28. The broken line describes the chemical freeze-out conditions of fixed total 
density of baryons plus antibaryons, rij, + rig = 0.12/fm 3 from Ref. 1491 . The full line 
represents the condition of the fixed energ y/particle ~ 1.0 GeV from Fig. J22). The 
freeze-out points are as in Fig. 1271 . 



experimental data from the top AGS up to RHIC energy. However, in the 
energy range from SIS to AGS it slightly overestimates the freeze-out tem- 
perature for a given chemical potential. Consequently, e.g. the yield of the 
strange/non-strange particle ratios obtained at SIS turns out to be too 
large. The freeze-out conditions determined by the extensive thermody- 
namical observables are in addition very sensitive to the size and the model 
that describes repulsive interactions between hadronic constituents. 

The condition of fixed (E)/(N) ~ 1 GeV, is very insensitive to repulsive 
interactions. Independently on how the rep ulsive interactions are imple- 
mented, that is through a mean field potentiai^-^, a n effect ive hard cor J5BI 
or a thermodynamically consistent implementation 1 ^ 1 8^ 1 the freeze-out 
line in Fig. (|27[) is hardly modified. However, the energy per particle is be- 
ing sensitive to the composition of the collision fireball. Considering heavy 
fragments like e.g. the He or Li as being the constituents of a thermal fire- 
ball would change the line shown in Fig. J2ZJ- In general these fragments 
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would make the line steeper below the AGS energy. This is an open question 
whether at the SIS energy, the multiplicity and the spectra of such com- 
posite objects are of thermal origin and can be reproduced with the same 
parameters as all other hadrons. For higher energies beyond AGS, however, 
the above is not excluded as discussed in Section 3.04 and in Ref. I|175fl . 

For the phenomenological determination of freeze-out parameters for 
different collision energies we use in the following the requirement (E) / (N) ~ 
1 GeV. 



5.1. Chemical freeze out and the QCD phase boundary 

The chemical freeze-out temperature, found from a thermal analysis 



of experimental data in Pb-Pb collisions at the SPS and in Au-Au colli- 
sions at RHIC energy is remarkably consistent, within errors, with the crit- 
ical temperature T c ~ 173 ± 8 MeV obtained-^ from lattice Monte-Carlo 
simulations of QCD at a vanishing net baryon density. Thus, the observed 
hadrons seem to be originating from a deconfined medium and the chemical 
composition of the system is m ost lik ely to be established during hadroniza- 
tion of the quark-gluon plasmaSEEl The observed coincidence of chemical 
and critical conditions in the QCD medium at the SPS and RHIC energy 
open the question if this property is also valid in heavy ion collisions at 
lower collision energies where the statistical order of the secondaries is phe- 
nomenologically well established. 

Recently, first attempts have been made to ext ent lattic e calculations 
into the region of finite [Lb ■ This provided an estimat(E3IZiof the location 
of the phase boundary at finite baryon density. The generic problem of the 
Monte-Carlo simulation of QCD with the finite chemical potential, related 
with a complex structure of the fermionic determinant, was partly over- 
come. The reweighting method, in which the physical observables at finite 
fiB are computed by simulating the theory at vanishing u h 1 21 wa s suc cess- 
fully applied and first results on the phase boundary were obtained^. The 
region of applicability of this approach and uncertainties on the results due 
to a small lattice size and large strange quark mass are still, however, not 
well established. Another efficient method, at least for low baryo n den sity, 
is based on the Taylor expansion in fis of any physical observabl e ! "^71 The 
coefficients of the series are calculated at vanishing \ib , and thus could be 
obtained using a standard Monte-Carlo method. This procedure was re- 
cently applied to get a series expansion of the critical temperature in terms 
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Fig. 29. A Comparison of the chemical freeze-out curve from Fig. 1271 with the phase 
boundary line. The upper thin line represents the LGT results obtained in Ref. 11771 
and the lower thin line describes the conditions of constant energy density that was fixed 
at fi = 0. The upper point with crossed error bars denotes the end-point of the crossover 
transition from Ref. 11781 . 



of the critical 

Fig. H29f) shows the results on the position of the phase boundary that 
were obtained using the methods indicated above together with the freeze- 
out curve from Fig. i|27|) . The upper thin- line represents an extrapolation 
of the leading, (hb)c order, ter m in the Taylor expansion of T c to a larger 
values of the chemical potentia l It is interesting to note that, within 
statistical uncertainties, the energy density along this line is almost con- 
stant and corresponds to e c ~ 0.6GeV/fm 3 (thin lower line in Fig. J2SJ0; 
that is the same value as found on the lattice at (1b = 0. It is thus, con- 
ceivable that the critical (fig — T c ) surface is determined by the condition 
of fixed energy density. This can be also argued phenomenologically. The 
transition from a confined to deconfined phase could appear if the particle 
(like in percolation models) or energy density is so large that hadrons start 
to overlap. It should not be important if this density is achieved by heat- 
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ing or compressing the nuclear matter. Thus, since the percolation type 
argument! ' ^ is well describing critical conditions at (1b = it could be 
also valid at finite [Ib- 

Fig. I|29|l shows that the chemical freeze-out points at SPS and RHIC 
energy are indeed lying on the phase boundary. The results of SPS at 40 
A-GeV and top AGS are already below the boundary line. However, it is not 
excluded that also at these lower energies the collision fireball in the initial 
state appears in the deconfined phase. The initial energy density expected 
at AGS is of the order of 1 GeV/fm 3 (see Section 1.1) thus, it is larger than 
the critical energy density along the boundary line in Fig. I|29ll . 

The canonical suppression effects for strangeness production, were shown 
to be negligible already at the top AGS energies. Here, strangeness was un- 
corrected and well described by the GC approach. It is quite possible that 
the asymptotic GC formulation and the maximal thermal phase space for 
strangeness is achieved if in the initial state the system was created in a 
deconfined, QGP phase. The abundant production of strangeness in the 
QQfEB together with a long range correlations during a non-perturbative 
hadronization results in strangeness population that maximizes the entropy 
in the GC limit. In this context the energy range between SIS and 40 A- GeV 
is of particular interest and it is expected to be covered by the plannecP^I 
for the future heavy ion experiments at GSI. 



6. Particle yields and their energy dependence 

The hadronic composition in the final state obtained in heavy ion colli- 
sions is determined solely by an energy per hadron to be approximately 1 
GeV per hadron in the rest frame of the system under consideration. This 
phenomenological freeze-out condition provides the relation between the 
temperature and the chemical potential at all collision energies. 

The above relation together with only one measured particle ratio, e.g. 
the ratio of pion/participant° as shown in Fig. (|30J) establishe d the 
energy dependence of the two thermal parameters T and /is . Consequently, 
predictions of particle excitation functions can be given in terms of the 
canonical statistical model. An alternative approach would be to interpolate 
and/or parameterize the energy dependence of the [1b and then using the 



°The mean number of pion multiplicity is defined as: (7r) = 1.5((7r + ) + {it }) whereas 
the number of participant is calculated as the number of wounded nuclcons 
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Fig. 30. The total number of pions per wounded nucleon ((7r)/A r w ) versus the center- 
of-mass energy. The data at lower energies in A— A as well as in p-p collisions are from 
Refs. fnmri55l . The RHIC results are from Ref. ITTSt . The short-dashed and dashed 
lines are a fit to the data. 



unified freeze-out condition of (E) / (N) ~ 1 to get the energy dependence 
of T. The energy dependence of the chemical potential was showrpSto be 
well parameterized as 

where a ~ 1.27 GeV and b ~ 4.3 GeV. The result of this parameterization 
is shown by the full line in Fig. together with the energy dependence 
of the freeze-out temperature. 

In the statistical approach, the knowledge of T(y/s) and /ib(V^) deter- 
mines the energy dependence of different observables. Of particular inter- 
est are the ratios of strange to non-strange particle multiplicities as well 
as the relative st range ness content of the system that is measured by the 
Wroblewski factoiPll. 

We turn our attention first to a study of the energy dependence of the 
Wroblewski ratio defined as 

2(ss) , , 

A -HTM' (n9) 
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Fig. 31. Behavior of the freeze-out baryon chemical potential hb (upper curve) and 
the temperature T (lower curve) as a function of energy from Ref . 1531 . The temperature 
T as a function of beam energy is determined from the unified freeze-out conditions of 
fixed energy/particle. 



where the quantities in angular brackets refer to the number of newly 
formed quark-antiquark pairs, i.e., it excludes all the quarks that are present 
in the target and projectile. 

The quark content used in this ratio is determined at the moment of 
chemical freeze-out, i.e. from hadrons and especially, hadronic resonances, 
before they decay. This ratio is thus not easily measurable unless one can 
reconstruct all resonances from the final-state particles. 

The results are shown in Fig. 1)32(1 as a function of center of mass en- 
ergy y/s. The values calculated from the experimental data at chemical 
freeze-out in central A- A collisions have been taken from reference l|4l)fl . p 
All values of A s were extracted from fully integrated data besides RHIC 
where the STAR collaboration results on particle ratios measured -*-^ at 
mid-pseudorapidity were used. The solid line in Fig. I(32JI describes the 
statistical model calculation^^ in complete equilibrium along the unified 



p Here the statistical model was fitted with an extra parameter 7 S to account for possible 
chemical undersaturation of strangeness. At the SPS, f s ~ 0.75 was required to get the 
best agreement with An data. See the discussion in chapter 2 concerning this issue. 
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Thermal Model with y =1 and radius = 7 fm 

1.1 - Thermal Model with II ^ = 0.0 - 

~ Thermal Model with |l = 0.0 and radius = 1.2 fm ~ 



k =2<ss>/<uu + d3> - 




Vs (GeV) 

Fig. 32. The Wroblewski ratio A s as a 
function of y/s. For the description of the 
lines see text. The points are the statistical 
model results calculated with thermal pa- 
rameters obtained from the fit to measured 
particle yields. 
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Fig. 33. Contributions to the 
Wroblewski factor from strange baryons, 
strange mesons and hidden strange 
particles. Full line is a sum of all these 
contributions. 



freeze-out curve and with the energy dependent thermal parameters shown 
in Fig. I|3ip, From Fig. l.'il'l) one sees that around 30 A- GeV lab energy the 
relative strangeness content in heavy ion collisions reaches a clear and well 
pronounced maximum. The Wroblewski factor decreases towards higher in- 
cident energies and reaches a limiting value of about 0.43. 

The appearance of the maximum can be related to the specific depen- 
dence of hb on the beam energy. In Fig. (|32[) we also show A s calculated 
under the assumption that only the temperature varies with collision energy 
but the baryon chemical potential is kept fixed at zero (dotted line). In this 
case the Wroblewski factor is indeed seen to be a smooth function of energy. 
The assumption of vanishing net baryon density is close to the prevailing 
situation in e.g. p-p and e + -e~ collisions. In Fig. (|32|l the resu lts fo r A s 
extracted from the data in p-p, p-p and e + — e~ are also included!^. The 
dashed line represents the results obtained with fis = and a canonical 




Fig. 34. The ratio of kaon to pion measured in heavy ion collision at different collisions 
energies. The left-hand figure describes midrapidity data, whereas the right hand figure 
represents the ratio of fully integrated yields. Data at SIS, AGS, SPS and RHIC energy 
are taken from Refs. (185 187 1901. The short-dashed line describes the statistical 
model predictions along the unified freeze— out curve. The right hand figure also shows 
the parameterization of the p— p data (full-lines) from Ref. 11911 and the canonical model 
results (dashed-line). 



radius of 1.2 fin. There are two important differences in the behavior of 
X s in elementary compared to heavy ion collisions. Firstly, the strangeness 
content is smaller by a factor of two. This is mainly because in the elemen- 
tary collisions particle multiplicities follow the values given by the canonical 
ensemble with radius 1.1-1.2 fm whereas in A- A collisions the grand canon- 
ical ensemble can be used, thus strangeness is uncorrelated and distributed 
in the whole fireball. Secondly, there is no evidence, at the moment, of a 
significant maximum in the behavior of A s in elementary collisions. 

The importance of finite net baryon density on the behavior of A s is 
demonstrated in Fig. 1)33(1 showing separately the contributions to (ss) com- 
ing from strange baryons, strange mesons and from hidden strangeness, i.e., 
from hadrons like <j> and rj. As can be seen in Fig. 1(33(1 . the origin of the 
maximum in the Wroblewski ratio can be traced to the contribution of 
strange baryons that is strongly enhanced in the energy range up to 30 
A-GeV. 

The appearance of the maximum in the strangeness content of the colli- 
sions fireball can be also justified on the level of different ratios that includes 
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strange particles. The measured ^^^ K + /n + ratio (see also Fig. (|34|l ) is 
a very abruptly increasing function of the collision energy between SIS up 
to the top AGS energy. At higher energies it reaches a broad maximum be- 
tween 20 A-GeV - 40 A-GeV and gradually decreases up to RHIC energy. In 
microscopic transport model J " ' the increase of the kaon yield with colli- 
sion energy is qualitatively expected as being due to a change in the produc- 
tion mechanism from associated to direct kaon emission. However, hadronic 
cascade transport models do not, until now, provide quantitative explana- 
tion of the experimental data in the whole energy range. This is evident 
in Fig. l|35jl where the comparisons of RQMD, URQMD and BUU micro- 
scopic transport models with experimental data is presented. The RQMD 
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Fig. 35. The ratio of kaon to pion measured in heavy ion collision at different col- 
lisions energies in comparison with microscopic transport models. The left-hand figure 
represents the RQMeEEH UrQMrJSHSSl and the statistical model predictions. The 
right hand figure shows the BUU results 192 . For the description of data see Refs. 1571 



provides a good description of the high energy data. The URQMD works 
quite well at the low energy, however, underestimates the yields at the SPS 
and RHIC energy. The BUU on the other hand is well suitable for SIS en- 
ergy range, however, shows different energy dependence than that obtained 
in experiments. The statistical model in the canonical formulatiorl^,on 
the other hand, provides a good description of the Kjix midrapidity ra- 
tio from SIS up to AGS as seen in Fig. 134(1 . The abrupt increase from 
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SIS to AGS and the broad maximum of this ratio are consequences of the 
specific dependence of thermal parameters on the collision energy and the 
canonical strange ness supp ression at SIS. A drop in the K + /ir + ratio for 
4tt yields reportecffEZUJMH ^ y the NA49 Collaboration at 158 A-GeV (see 
Fig. (|35|l is, however, not reproduced by the statistical model without fur- 
ther modifications, e.g. by introducing an additional parameter j s ~ 0.7^^ 
that accounts for additional suppression of strang eness . We also note that 
an abrupt drop in the K + /n + ratio is predictecP^ in a special model 
with particular conditions on the early stage of the collisions. This model, 
however, neglects completely the production of strange particles from the 
hadronization of gluons. 




Fig. 36. Particle ratios in A-A collisions versus the center-of-mass energy. Data at 
the SPS are fully integrated NA49 results. Th e cor responding ratio at the top AGS 
was obtained from E810 results on H~ measured in Si— Pb collisions in the rapidity 
interval 1.4 < y < 2.9, normalized to the full phase-space values of 7r+ and K~ yield 
obtaineoEElnSi— Au collisions by E802 Collaboration. The lines represent statistical 
model resultJ^SI along the unified freeze-out curve. 



The appearance of the maximum in the relative strange /non-strange 
particle multiplicity ratios already seen in K + /n + is even more pronounced^ 
for strange baryon/meson ratios. Fig. Ij^fijl shows the energy dependence of 
A/7r + and S~/-7r + . There is a very clearly pronounced maximum, especially 
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in the A/ir + ratio. This maximum is related with a rather strong decrease 
of the chemical potential coupled with only a moderate increase in the as- 
sociated temperature as the energy increase*^. The relative enhancement 
of A is stronger than that of H — . There is also a shift of the maximum to 
higher energies for particles with increasing number of strange quarks. This 
is because an enhanced strangeness content of the baryon suppresses the 
dependence of the corresponding ratio on fig. This is also seen for S~ /K~ 
ratio that shows a substantially narrower maximum since the strangeness 
dilution effect is compensated by the strangeness content of the K~ . The ac- 
tual experimental data for both A/tt + and "Er /tt + ratios shown in Fig. (|36|) 
are following the predictions of the statistical model. However, as in the case 
of kaons, midrapidity results are better reproduced by the model than 4tt 
data. 
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Fig. 37. The left-hand figure, shows the ratio of K + /K~ as a function of (ir + — 
n°) /Apart - Points are the experimental results, a line is the statistical model result along 
the unified freeze-out curve. The right-hand figure shows K~ /K+ ratio that appears to 
be constant as a function of centrality from SIS up to RHIC energy. Data are from the 
STAR, NA49, E866, and KaoS Collaborations. The broken lines are statistical model 
results. 



The statistical model predicts that if at least two different ratios of 
non strange particles are constant with centrality then also strange par- 
ticle/antiparticle ratio should be centrality independent. Dynamically this 
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is a rather surprising result as strange particles and their antiparticles are 
generally produced and absorbed in the surrounding nuclear medium in a 
different way. This is particularly the case at lower energies (up to AGS) 
where e.g. K + and K~ are predominantly produced due to ttN — > AK + 
and it A — > K~ N processes. In addition K + , mesons feel a repulsive poten- 
tial whereas K~ mesons are attracted. Thus, the prediction of the thermal 
model that the K + /K~ ratio is centrality independent was dynamically 
unexpected. Figure (|3"?|l represents the energy and centrality dependence 
of the K + /K~ ratio from SIS to RHIC energy. The statistical model pre- 
dictions are seen in Fig. H37l -right) to agree remarkably well with the data. 
The results of Fig. (|37l -right') could be considered as the evidence of an 
apparent chemical equilibrium population of kaons in the final state. This 
behavior of data is also seen on a different level. The chemical equilibra- 
tion of the associated production of K + with a hyperon and strangeness 
exchange production of K~ , indicated ab ove, should result in a linear de- 
pendence of the K + /K~ ratio on i:/A par l ^ \ Fig. 1)371 - left) shows that, 
indeed, in the energy range from SIS up to AGS, and almost independently 
from the colliding system, the above prediction is valid. It is also clear from 
this figure that between AGS and SPS the production mechanism of strange 
mesons is changing. 

The results for the K + and K~ excitation function show an interesting 
behavior when expressed as a function of available energy t/s — y/sth (see 
Fig. El-left j!3 

The threshold energy y/sth corresponds to the production 
threshold in N-N collisions. For K + mesons ^Jsth — 2.548 GeV, whereas 
for K~ the corresponding value is ^Jsth = 2.87 GeV. It turned out that in 
this representation the measured yields of K + and K~ , close to threshold, 
in heavy ion collisi ons a re about equal while they differ in p-p collisions by 
a factor of 10 - 10Cp22l. Fig. JJHl-left) shows the experimental data together 
with the canonical statistical model resultg32l. In the range where \fs — ^/s t h 
is less than zero, the excitation functions for K + and K~ , obtained in the 
model, cross each other, leading to the observed equality of K + and K~ 
at SIS energies. The yields differ at AGS energies by a factor of five. The 
difference in the rise of the two excitation functions can be understood 
in the statistical model as being due to the different dependence of K + 
and K~ yields on fis- The K + yield is strongly fis dependent through 
associated production with A whereas K~ yield is not directly effected by 
[iB- Consequently, the excitation functions, i.e. the variation with T, exhibit 
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a different rise for kaons and anti-kaons. 




Vs-Vs^ (GeV) 

Fig. 38. The left-hand figure: calculated /A par t and K~ /A par t ratios of 
yield/participant in the statistical model as a function of ^ /s — ^/sth f° r Ni+Ni col- 
lisions. The points are the res ults f or Ni+Ni collisions at Slj 119 ^ 2U l and for Au+Au 
collisions at 10.2 A- GeV (AGSjI^U energies. T he rig ht-hand figure: the yield of 4>/K 
ratio in Ni-Ni and Au— Au collisions calculated I in the statistical model along the 
unified freeze— out curve. The data points are from Ref. ( 203 1 

In the statistical model all particle species are considered to be on shell. 
Thus, the statistical model reproduces the kaon yields and their excitation 
functions with current particle masses. A transport calculation^^E^EHl 
on the other hand, required in-medium modifications (a reduction) of the 
K ~ mass (as expected for kaon in the nuclear medium) in order to describe 
the measured yields. These differences are, however not necessarily in con- 
tradiction, as the transport model describes a time evolution of particle 
production, whereas the statistical models are only valid for particles that 
are measured in the final state, where the on-shcll conditions are to be ex- 
pected. 01 It is conceivable that the apparent chemical equilibration observed 
in the data at SIS energies could be a direct consequence of in-medium ef- 
fects. One such possibility is an increase of the in-medium production cross 



q Nevertheless, in transport models the kaons and anti-kaons are produced in different 
time during the collision, thus also at different temperatures. In the statistical model 
there is a common temperature for kaons and their antiparticlcs. 
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section of strange particles. 

In the context of particle production in heavy ion collisions a particular 
role has been attributed to the vector meson resonances. The measurements 
of these particles could possibly provide an information on the chiral sym- 
metry rest orati on ' ' ' ' ' and in medium effects due to the collision 
broadening ]!^ ! of their decay widths. The production of cf> meson is of par- 
ticular interest due to its ss content that could make it a very sensitive 
probe of strangeness production in the early stage. The yield of ^ mesons 
at the SPS was found to be very different when it was reconstructecP^EHI 
from the K + K~ or /i + fi~ deca y ch annels . Possible scenarios of this differ- 
ence were extensively discusse d in the literature. In view of the above 
problems it is interesting to check if the production of 4> mesons follows the 
statistical order. Figure 1381 -right) sh ows e xpe rimen tal values for the ratio 
of (f)/K- obtained at SIS 203 206 , AG^OHl, SP^Hland RHIC energies^. 
These results are compare d 2 ^ with the canonical model for Ni-Ni and 
Au-Au collisions. The 4>/K~ ratio is seen in Fig. 1)381 -right) to exhibit a 
broad maximum as a function of collision energy that appears around 40 
A-GeV. Thus, the maximum strangeness content indicated in Fig. I|36(l is 
also visible here. It is surprising that the yield of <fi/K~ within 40% is the 
same in a very broad energy range from AGS up to RHIC. The Ni-Ni re- 
sults seems to follow this systematics. Howe ver, we have to point out that 
the most recent, still preliminary result f the FOPI Collaboration in- 
dicate a much larger value of the ratio than it is shown in Fig. J3H1 -right). 
Also at the SPS the resultJ223l f the NA50 Collaboration, extrapolated 
with the measured slope to lower p t , would give the value that is by a 
factor of 3 - 5 larger then that obtained^] by NA49 Collaboration. The 
above uncertainties notwithstanding, the statistical model follows the trend 
of data and reproduces their magnitude remarkably well. It is interesting to 
note the differences in the canonical model on the energy dependence of the 
<j)fK~ for Ni-Ni and Au-Au collisions. The Ni-Ni line in Fig. -right ) 
shows a sharp minimum at 3 < ^fs < 4 that appears due to a strong de- 
crease of the K~ yield through the canonical suppression. In the canonical 
model, applied here 1 the <f> meson yield is unconstraint. This would not be 
the case if the strangeness undersaturation factor, instead of the correlation 
volume, would be implemented. 



r The results were obtained in the canonical model without strangeness suppression factor. 
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Another interesting future of the 6 yi eld i s its centrality dependence. 
The ratio of (j)/K~ measured at AG^DDI SP^HIl and RHICpHl energies 
is centrality independent whereas both <j> and K~ are showing a strong 
centrality dependence that increases with decreasing beam energy. Thus, 
the 4> meson shows a similar centrality dependence as kaons. 

The centrality dependence of the K/tt ratio at the SPS was recently 
described I I, in the context of a statistical model, as being due to cen- 
trality dependence of strangeness undersaturation parameter A s . This is 
quite a conceivable scenario. The K + /K~ as well as p/p ratios are only 
very weakly centrality dependent, thus the temperature and baryon chem- 
ical potential, should be rather constant with A part . The change of the 
K/t: ~ X s with A part could be thus attributed to A s . The same arguments 
could be also applied to the AGS and RHIC results. However, the 4>/K~ is 
to be proportional to A s , consequently the <f>/K~ ratio should show similar 
Apart dependence as K/ir ratio. The above, however, is not confirmed ex- 
perimentally. This discrepancy calls into question the concept of a 7 S factor 
as the relevant parameter characterizing strangeness production in heavy 
ion collisions (see also the discussion in Section 2). 



7. Lifting of the strangeness suppression in heavy ion 
collisions 

An enhanced production of strange particles compared to the suppressed 
strangeness yi eld obser ved in collisions between elementary particles was 
long sugg as a possible signal of the QGP formation in heavy 

ion collisions. 8 In the QGP the production and equilibration of strangeness 
is very efficient due to a large gluon density and a low energy threshold 
for dominant QCD processes of ss productior l^ l ^ l In hadronic systems 
the higher threshold for strangeness production was argued in Ref. to 
make the strangeness yield considerably smaller and the equilibration time 
much longer.' 



s Strangeness enhancement in heavy ion collisions was recently discussed and interpreted 
at the parton level in Ref. {211 1 without requirements of the QGP formation in the initial 
state. 

'Recently, it was argued in Rcfs. 1 128 33 1, that multi-mesonic reactions could accelerate 
the equilibration time of strange antibaryons especially when the hadronic system is hot 
and very dense. This argument does not apply, however, for strange baryon, where also 
strong enhancements are seen. 
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Based on such arguments predictions have been developed for experi- 
mental signatures of deconfmement. Key predictions ar d ' \ 

i) the disappearance of the strangeness suppression observed in collisions 
among elementary particles leading to e.g. an enhancement of multistrange 
baryons and anti-baryons in central A-A collisions, with respect to proton 
induced reactions. 

ii) chemical equilibration of secondaries: the appearance of the QGP be- 
ing close to a chemical equilibrium and subsequent phase transition should 
in general drive hadronic constituents produced from hadronizing QGP to- 
wards chemical equilibrium. _____ 

Heavy ion experiments at the CERN SPS report ec r 1 1 '4 \ actually a 
global increase of strangeness production from p-p, p-A to A-A collisions. 
This effect is seen e.g. in Fig. (|__ -right) that shows the enhancement of 
the K + /tt + ratio in Pb-Pb relative to p+p collisions. There is indeed an 
increase by a factor of about two in strangeness content when going from 
p+p to heavy ion collisions. This is also seen on the level of the Wroblewski 
factor in Fig. (J_J. 
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Fig. 39. Particle yields per participant in Pb- Pb r elative to p— Be and p— Pb collisions 
centrality dependence. The data are from WA9 ~ 3 and NA5 ~ _ Collaborations. 

A large strangeness content of the QGP plasma should, following (ii), 
be reflected in a very specific hierarchy of multistrange baryon ~ 5 implying 
an enhancement of Vl > 3 > A. Fig. shows the yield/participant in 
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Pb-Pb relative to p-Be and p-Pb collisions measure the WA97 

and NA57 Collaborations. Indeed the enhancement pattern of the antihy- 
peron yields is seen to increase with strangeness content of the particle and 
in WA97 data there is a saturation of this enhancement for N woun d > 100. 
The recent results of the NA57 collaboration are showing in addition an 
abrupt change of the anti-cascade enhancement for a lower centralit y. Sir n- 
ilar behavior was previously seen on the level of K + yield measured I at 
Slab = by the NA52 experiment in Pb-Pb collisions. These results are very 
interesting as they might be interpreted as an indication of the onset of a 
new dynamics. However, a more detailed experimental study and theoreti- 
cal understanding are still required here. Until now there is no quantitative 
understanding of this exceptional threshold behavior of EE. 

A number of different mechanisms were considered to describe the mag- 
nitude of the enhance ment and cent rality dependence of (multi)strange 
baryons measure d^ 1 $9 1 20 . 9 1 2 1 5 I-ED1 by the WA97 Collaboration. Stud- 
ies using microscopic transport models make it clear that data shown in 
Fig. (|39J) can not be explained by pure final state hadronic interactions. 
The combination of the former with additional pre-hadronic mechanisms 
like sting formation and their subsequent hadronization, baryon junction 
mechanism, color ropes or a color flux tubes overlap improves the agree- 
ment with the measured enhancement pattern and the magnitude for the 
most central collisions. However, the detailed centrality dependence is still 
not well reproduced within the microscopic models. 

The results of Section 2 have shown that the statistical model provides a 
satisfactory description of strange and multistrange particle yields in A-A 
collisions. The midrapidity data were argued in Section. 1 to be reproduced 
by the statistical model in a full equilibrium. The fully integrated results, on 
the other hand can be successfully reproduced when including a strangeness 
undersaturation parameter j s ~ 0.75 that accounts, at the SPS, for a 25% 
deviation from equilibrium. 

Strangeness production in p-p collisions was shown in Section 5 to be 
consistent with the canonical statistical model. The abundance of single- 
strange particles could be described by this model by including, as in heavy 
ion collisions, the strangeness undersaturation factor 7 S ~ 0.51 or a corre- 
lation volume V ~ V p that accounts for the locality of strangeness con- 
servation. Consequently, the enhancement from p-p to A-A collisions of 
strangeness 1 particles could be well described in terms of the statistical 
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model as the transition from the canonical to the asymptotic GC limi1@Sl. 




Fig. 40. The left hand figure: the canonical suppression factor F s ~ I s (Vx)/ Iq(Vx) 
as a function of V ~ (N pa rt) with the argument x calculated for a hadron resonance gas 
at T = 168 MeV and (j,b = 266 MeV. The right hand figure: statistical model results 
on centrality dependence of the relative enhancement of Air particle yields/participant 
in central Pb— Pb to p-p collisions expected at y/s ~ 17 GeV. 

One of the consequences of the canonical model is a very particular vol- 
ume dependence of multistrange particle densities. In heavy ion collisions 
this volumes is usually assumed to scale with the number of participants. 
From Eq. (|97|l it is clear that the canonical suppression should increase with 
strangeness content of the particles. Indeed the approximate strangeness 
suppression factor F s ~ I s (Vx) /Iq (Vx) for a fixed Vx is a decreasing func- 
tion of s. This is particularly evident in the limit of (Vx) << 1 where 
I s (Vx)/Iq(Vx) ~ (Vx) s . The suppression factor is quantified in Fig. 1(101 
-left) that indicates the expected suppression pattern. 

Fig. I|41l-left) shows predictions-*^* of the canonical model for the multi- 
plicity/participant of f2, S, and A relative to their value in p-p collisions. 
Thermal parameters, T = 168 MeV and hb = 266 MeV, that are appro- 
priatd^ for a description of central Pb-Pb collisions at 160 A-GeV, were 
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used here and assumed to be centrality-independent." The canonical vol- 
ume parameter was taken to be proportional to the number of projectile 
participants, V ~ V N palt /2 where V = 4nR^/3 with R p = 1.15 ± 0.5. 
The volume Vq is then, depending on the particular choice of R p , of the 
order of the volume of a nucleon. Figure (|40l -right) indicates that the 
statistical model in the canonical ensemble reproduces the basic features of 
WA97 data: the enhancement pattern and enhancement saturation for large 
^4part ■ The appearance of the saturation of the enhancement indicates that 
the grand canonical limit was reached. It is also clear from Fig. J2]-left) 
that this saturation is shifted towards a larger centrality with increasing 
strangeness content of the particle. 



1000 



: oo 



1 




















a 

*e 
















> s 










• E9 7 7 




























■ ES66 


z / 




\ 








♦ NA49 


/' 














/ 
i 




*e 




\ 






i 




\ 




\ 




















i 














! / 








< 






if 




























a 




2 










a/ 














if' 






























6 8 


10 


12 14 16 18 



<N„ 



100 
„> 



1 000 



[GeV] 



Fig. 41. The left hand figure: the statistical model results on the centrality dependence 
of the relative enhancement of A, S and Q yields/participant in central Pb— Pb to p— p 
collisions at y/s = 8.73 GeV. The right hand figure: K + /n + ratio in A+A relative to p+p 
collisions. For the compilation of data, see Refs. (191 190). The dashed line represents 
the statistical model results. 



The essential predictiorPai of the canonical statistical model is that the 
strangeness enhancement from p-p to A-A collisions should increase with 
decreasing collision energy. This is a direct consequence of e.g. Eq. (|§7|) 



"For a detailed discussion of the possible centrality dependence of thermal parameters 
sec e.g. Ref. 12091 . 
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where the canonical suppression factor is seen to be a decreasing func- 
tion of temperature and thus also collision energy. Fig. (|4 1 1 -right) shows 
the compilation of data on the K + /ir + ratio in A-A relative to p-p col- 
lisions from Ref. I|191l [190 ) . This double ratio could be referred to as a 
strangeness enhancement factor. The enhancement is seen in data to be 
the largest at the smallest beam energy and is decreasing towards higher 
energy. The line is a smooth interpolation between the canonical model 
results for y/s = 17.3,12.3,8.73,5.56 GeV, calculated with the thermal 
parameters T and /i^ that were extracted from Fig. I|27|) . The canonical 
volume parameter was taken the same as used in Fig. (|16t right). The en- 
hancement seen in Fig. (^2 -r ight) is due to the suppression of the K + /ir + 
ratio in p-p collisions with decreasing energy and not due to a dilution of 
this ratio by excess pions in the A-A system. The K + /n + ratio is known 
experimentally not to vary within 30% in th e energy ran ge from -y/s ~ 5 
GeV at AGS up to Js = 130 GeV at RHIcPZHWHHl 

Recent data of 

NA4d333 an d CERES^HEIIBIIl 

Collaborations on A yields exhibit a similar 
energy dependence of the enhancement factor as is seen in Fig. Q411 right) 
for kaons. 

The canonical model (see Eq. JEJl) also predicts that the multistrange 
baryon enhancement from p-p to A-A should be larger at lower collision en- 
ergy. Fig. Q411 left) shows the canonical model results for the strangeness en- 
hancement and its centrality dependence for Au-Au collisions at 40 A- GeV. 
The qualitative behavior of the enhancement is like that at the SPS. How- 
ever, the strength of the enhancement is seen to be substantially larger. For 
VI it can be as large as a factor of 100. The enhancement of S in A-A relative 
to p-p can be deduced from data. Indeed, combining the Si-Pb results for 
S~ production obtainea^^^ by E810 Collaboration and the Si-Au results 
for pion or K~ yields obtainecP^ by E802 Collaboration in collisions at 
the top AGS energy, one can estimate that E~ /tt + ~ 0.0076. Within errors 
this agrees with the value of E~ /tt + ~ 0.0074 obtainecPSEI 

by NA49 in 

Pb-Pb at yfs = 17.3 GeV, as seen in Fig. I|36|) . In p-p collisions the S~/7r + 
ratio is obviously a strongly decreasing function of beam energy. From the 
above one would therefore expect that the relative enhancement of S~ from 
p-p to A-A collisions should be larger at AGS than at SPS energy. This 
is seen in Fig. (|41l -leftV In addition ratios containing only newly produced 
particles, such as e.g. S~ /K~, are also larger at AGS than at SPS energy, 
see Fig. I|36[) . The above heuristic argument should be, however, verified by 
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more complete experimental data. 

Recently, the production of S~ was in the relativistic trans- 

port model in the energy range from 2-11 A-GeV . The authors discussed 
the equilibration and the influence of a phase transition on 3 yields in A-A 
collisions. There, it was argued that, for beam energies above 4 A-GeV, 
there should be a sharp increase of the S~ yield if there is a deconfinement 
transition in the collisions fireball. This is an interesting conjecture. How- 
ever, with the presently available experimental data this prediction cannot 
be tested. The canonical statistical model, predicts continuous increase of 
3 yield with collisions energy as it is seen in Fig. i|36|) . 

The prediction of the statistical model on the energy dependen ce of 
strange particle yields is also in contrast with the UrQMD resultJ^^. 
There, the production of strangeness is very sensitive to the initial condi- 
tions. In UrQMD the early stage multiple scattering may imply an inc rease 
of the colour electric field strength due to an overlap of produced string 
Consequently, according to the Schwinger mechanism, this should increase 
the production of (multi)strange baryons. Under similar kinematical con- 
ditions as at the SPS, the UrQMD model predictJ2HI at RHIC an increase 
in relative strength of £1 yields from p-p to A-A by a factor of 5. A recent 
analysis of multistrange baryon yields in A u-A u collisions at RHIC energy 
in the context of Dual Parton Model, lead J^^ to a smaller increase of the 
enhancement at RHIC energies. 

The interpretation that the strangeness enhancement is explained by 
canonical effects is not necessarily in contradiction with the heuristic predi- 
ctionJSl that strangeness enhancement and its pattern are due to a quark- 
gluon plasma formation in A-A collisions. This is particularly the case 
if one connects the asymptotic grand canonical description of strangeness 
production in A-A collisions with the formation of a QGP in the initial 
state. 



8. Conclusions and outlook 

This article has focussed almost exclusively on the use of the statistical 
approach to understand yields of different particle species that have been 
measured in heavy ion collisions. We have discussed a statistical description 
of the conservation laws and described their kinetic implementations. We 
have argued, both on the qualitative and quantitative level, that exact con- 
servation of quantum numbers is of crucial importance when applying the 
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methods of statistical physics in the context of heavy ion and particle colli- 
sions. We have presented the systematics of particle production in heavy ion 
collisions from SIS up to LHC energy and discussed particular properties of 
strangeness production. One of the most intriguing results that comes from 
these investigations is the observation that particles seem to be produced 
according to the principle of the maximal entropy. In a very broad energy 
range, from y/s ~ 2.5 up to 200 GeV per nucleon pair, hadronic yields and 
their ratios observed in heavy ion collisions resemble a chemical equilibrium 
population along a unified freeze-out curve determined by the conditions of 
fixed energy/hadron ~ lGeV or complementary above SIS energy by fixed 
total density of baryons. Strangeness production follows this systematics 
from low to very high energy. However, there are some characteristic fea- 
tures of the system at chemical freeze-out in high energy central A+A col- 
lisions regarding strangeness production that are not present in low energy 
heavy ion and collisions among elementary particles. In nucleus-nucleus 
collisions, strangeness is un-correlated and redistributed in the macroscopic 
volume of a collision fireball and is conserved on the average. In hadron- 
hadron collisions the thermal phase space available for strange particles 
is strongly suppressed since, with only few particles produced per event, 
strangeness is strongly correlated in a volume that approximately coincides 
with the size of a nucleon, i.e. a distance over which color is confined. Thus, 
following the statistical kinetics, strangeness has to be conserved exactly 
and locally. The associated production and locality of strangeness conserva- 
tion is, in the context of the statistical model, the origin of the suppression 
of the thermal phase space for produced particles. Within this context the 
strangeness suppression observed in collisions among elementary particles 
finds its natural explanation. We also note that the suppression increases 
with the strangeness content of the particle as well as, for A-A collisions, 
with decreasing collision energy. A further consequence of the transition 
from the canonical to the grand-canonical regime is that strangeness pro- 
duction should be enhanced in A-A collisions compared to p-p collisions. 

At SPS and RHIC energies the freeze-out points approach the calcu- 
lated QCD phase boundary. This fact lends strong support to the inter- 
pretation that the matter produced in nuclear collisions at SPS and RHIC 
energies was first thermalized in the deconfined quark-gluon plasma phase 
and subsequently expanded through the phase boundary into a thermal gas 
of mostly elastic and quasi-elastic interacting hadrons. The above connec- 
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tion between the QCD phase boundary and the observed chemical freeze- 
out points is sometimes called into questiorP^l since hadron production 
in e + e~ and p-p or p-p can also be described in thermal models yielding 
a (apparently universal) temperature T e w 170 MeV. While the fact that 
T e is close to T values determined for heavy ion collisions at top SPS and 
RHIC energies might indeed reflect the fundamental hadronization scale 
of QCD, we have already noted above that there is an essential difference 
between thermal descriptions of central heavy ion collisions and elemen- 
tary particle reactions. Strange particle densities and their ratios can be, 
for heavy ion collisions at full AGS energy and higher, well described in the 
grand-canonical ensemble. In contrast, for a description of particle produc- 
tion in elementary collisions, local quantum number conservation, that is 
the canonical description is need. Consequently, in central nucleus-nucleus 
collisions at ultra-relativistic energies, strangeness percolates freely over 
volumes of thousands of fm 3 ! whereas in elementary processes is approxi- 
mately restricted to the size of nucleon. At top SPS and RHIC energies it 
is natural to conclude that in nucleus-nucleus collisions the percolation has 
its origin in the quark-gluon phase, lending further strong support to the 
interpretation that the "coincidence" between experimentally determined 
chemical freeze-out points and the calculated phase boundary implies that 
a deconfined phase was produced in such collisions. 

In the context of statistical physics the fact that the measured particle 
yields coincides with a thermal multiplicities calculated with a given statis- 
tical operator is a necessary and sufficient conditions to maintained ther- 
malization of the collisions fireball. In the sense of Gibbs interpretation of 
thermodynamics this implies that the parameters T and fiB reflect the ther- 
mal properties of the fireball. Furthermore in heavy ion collisions there are 
experimental observables indicating the appearance of thermodynamical 
pressure and correlations that are expected in a thermalized medium. The 
build-up of pressure and collectivity is seen in heavy ion collisions on the 
level of particle transverse momentum distributions and elliptic flow param- 
eters. An increase of average particle transverse momentum with particle 
mass, seen from SIS up to RHIC energy, is a typical property of a transver- 
sally expanding thermal medium. The appearance of strong elliptic flow is 
an indication of thermodynamical pressure that develops in the early stage 
in the collision. Finally, the measured particle number fluctuations are con- 
sistent with thermal expectations. Taken together, these observations lend 
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strong support to the thermodynamic interpretation of T and \ib with the 
concomitant QGP interpretation. 
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